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ABSTRACT 


An  extensive  investigation  is  presented  into  the  applicabil¬ 
ity  and  accuracy  of  the  method  of  aerodynamic  influence  coeffi- 
c,.ents  for  calculating  the  airload  distribution  on  thin  wings  of 
arbitrary  planform,  executing  small  steady  or  simple  harmonic 

ln  *  »upeyaonic  flow  of  perfect  gas.  Several  represen- 
4eXain5leS‘  invo}vin§  wings  with  various  combinations  of 
subsonic  and  supersonic  edges,  are  worked  out.  Where  possible. 

",ade  with_results  of  more  exact  linearized  theory. 
ThK^fe}^ive  *dva"ta8es  of  three  types  of  elementary  area  for 
subdividing  the  wing  planform  are  studied.  The  type  which  ap¬ 
pears  most  satisfactory  from  considerations  of  versatility  ac¬ 
curacy  and  simplicity  is  the  "Mach  box,"  a  rectangular  a^a 
element  with  diagonals  parallel  to  the  Mach  lines.  On  the  basis 
of  all  available  evidence,  the  method  is  concluded  to  be  satis - 
f°r  )J8e.in  flutte*  prediction  or  similar  applications, 

rtf-HC^n,nen^at^0r?  are  put  £orth  regarding  preparation  of  tables 
of  aerodynamic  influence  coefficients.  A  set  of  working  rules 
to  assist  the  engineer  in  using  these  coefficients  is  published 
in  a  subsequent  report  (Ref.  31).  K 
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State .(/?  =  <?  )  7  iii 
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LIST  OF  SYMBOLS 
Reference  serai-chord 
Streamwise  dimension  of  a  box 
Reference  length 

Amplitude  of  translational  oscillation 

v/7 


cho?rd(^/“)y  bi,ed  °n  referen« 

o?drbodxf7Tn»d  *tre*",'1,e 

Modified  reduced  frequency  based  on 

(  4  =  4 

Dimensionless  lift  per  unit  span 
Indices  designating  location  of  receiving  box 
Dimensionless  moment  per  unit  span 
Complex  amplitude  of  pressure  difference 

<A  *  AA  -  £»-£  J 

Dimensionless  amplitude  of  pressure  difference 
Pressure  contribution  from  row  ^7 

“‘SS)]  coordinates  [Eqs.  (2.33)  and 
Time 


Complex  amplitude  of  dovnwash  (positive  down) 
tribution  (8trength>  of  downwash  dis- 

Cartesian  coordinates 


Deflection  of  the  mean  surface  of  a  wine 
(positive  down)  6 


°Z  elementary  area  influencing  (* ,~) 
within  the  forward  Mach  lines  emanating  from 
'  / 


Necessary  functions  for  tabulatine 
square  boxes  (Eqs.  (2.9)] 


Necessary  functions  for  tabulating 
Mach  boxes  [Eqs.  (2.16)] 


AIC’s 

AIC's 


for 
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LIST  OF  SYMBOLS  (Con't.d) 

Aerodynamic  influence  coefficient  (Abbreviation) 
Pressure  influence  coefficients  associated  with 

£?"<Fl£r“h  di,trlb”tl#n  (*.8)  nd 

Aerodynamic  Influence  coefficient!  for  charac- 
teristic  boxes  [Steady  state,  Eqs.  (2.43)] 

-Complete  elliptic  integral  of  the  second  kind 
Fresnel  integrals  [Eq.  (2.43)] 

?!rnehr5preaentinS  fhe  Pressure  at  (  z,*  )  due 
to  unit  downvash  at  (  )  [Eq.  (6.1)] 

Functions  defined  by  Eqs.  (4.16) 

Bessel  function  of  the  first  kind  of  order /> 

Inverse  kerne i  representing  the  downvash  at 
'  X.fr  )  due  to  unit  pressure  at  ( )  [Eq.(6.2)] 

pe"J«l"le“  t°“1  Uft  *”d  to“l  “>««  re,- 

*?«C?l3C6?erlV‘tlVe*  (Cu,t0Mry  notation, 
Mach  number 

?EqCr*6  5)fS#me  influence  coefficient 

»Sfe(8t?Ur“  COefflcUn,:  fOT  square 

b«eJ3T^1C(|nf[yfnce  for  Mach 

Surface  of  planform 

Speed  of  flight  (Supersonic) 

Step  function  [Eq.  (6.4)] 


Angle  of  attack 

Amplitude  of  pitching  oscillation 

\/m^T 

Modified  reduced  frequency  based  on  the  side 
length  of  a  rhombus  [Eq.  (2.35)] 
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Indices  designating  location  of  sending  box 


Indices  designating  location  of 
relative  to  receiving  box 


sending  box 


Running  variables 
system 


in  tne  Cartesian  coordinate 


Density  of  the  undisturbed  medium 


Complex  amplitude  of  velocity  potential 
Circular  frequency  of  oscillation 
Modified  reduced  frequency  based  on  / 

(  a?  -  zk  nyp*) 

Sweep  angle  of  leading  edge 
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SECTION  I 
INTRODUCTION 

harmonicallyhoscillating^thin8  af™*?10?5™*1"10  forces  on  a 
in  supersonic  How  hJs  ?ece^ed  c^dl^M*  Uftlng  Surface 
the  original  work  of  von  Borbely  (Ref  aJtenti°n  since 

Two  distinct  approaches  have  h»l  f!n  1  j  ?  tbe  early  19^0’s. 
tions  for  various  wing  planforms.  °  l0Wed  in  develoPing  solu- 

Cions  for  eh«eilne«ised‘‘differeStU*CeqS«io!ldinfKe“Ct  solu‘ 

finite  wings  severff  Z  t*  aPPears  in  Refs.  2-5*  por 

ted  to  spefific  pleJfoS,  "iJne(?efrn!'»)SrlUtJ°3S  llml' 
quarter-infinite  wing  and  extend!*  r  3,4(  treated  the 
forms  with  oblique  straight  le£d?n»hfd  for"W*llo,\ t0  Pl«n- 
the  prescribed  motions  arearbitilfv dfhi*JRef'  5*  ’  Although 
the  velocity  potential  are  in  integral  p6  resu^ts  obtained  for 
tion  of  thii  theory  Is  lSHiSd  bJ  Slit&i?'  and«Jh*  aPPlic- 
cases  where  the  mode  of  oscillation  ?StfcJ1  considerations  to 
snanwise  coordinate.  Other  ron£i?k  f  independent  of  the 
6),  Stewartson  (Ref.  7)  and  Rott  C'8'i  Goodn,an  (Ref. 

solutions  for  deformations  wh?oJ  ill  8J .  drived  at  similar 
direction  only.  Ewlrd?"  eaSfvIi  SI  variable  in  the  stream 
Stewart  and  Li  for  the  oscillat-i™  *"****  concePt  was  used  by 
and  by  Chang  for  an  oscillliiS  XeI!C£j!f?Ur4win8  (Ref •  9) 
sonic  leading  and  trailing  edges  (Ref  ?ing  with  suPer" 

quently  found  that  thi.  cotJm?  ii  lfiil0> XJ  Wa?  subse- 
quencies.  An  extension  (Ref  i i \ 3  valid  only  for  low  fre > 

Li,  which  would  validate  the  method*?  pr£P°£ed  by  Stewart  and 
Concurrently  Watkins  treated  bhihn!L50r  h.igher  frequencies, 
in  rigid-body  motions  by  I  s«ies  «^f  J"  “I"*/  °«Ulatlng 
function  of  the  governing  integral  L,!!?  l0”,°f  che  k«nel 
frequency  (up  to  the  third  powl?  in  Ref Ci?o  in.terms  of  the 
seventh  power  in  Ref.  13).  P  4  Re^'  and  UP  to  the 

supersonic ^edges  performing  rigldrbodvUl*rt  °r  <ielt*  wing  '»ith 
Miles  (Ref.  if).  Proehlich8(Re!  in?  ?  7’S1?ns  wfs  studied  by 

r.l_i!ed_-for_ce_  expression,  for  this  "w?de°  SeU^essociStir*' 
aerodynamic^ theories  &rJun««£T4££u'fl!£* 

Te“lRepeoerred  ^  'he  1955  f0r  Publication  as  a  WADC 
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numerical  coXtatills  Vf  L°f  the  Piousness  of  the 
rigid  body  m^tions^nd'Ref^  aO^for0?^*?61^  *36  Ref*  19  for 
resentable  by  a  quadratic  form  y  harmonlc  deformation*  rep¬ 
low  frequency  and  thl  inlrffl  H°w*v*r»  the  limitation  to 

represent lng^nany  ?™es of  3e?nL^( th*  "u*dr*tlc  fon.  for 
ness  of  the*results^f  Ref.  M  ™  re,erl«  the  useful- 

for  [-odern^ircraft  w*™e£nd°t»il‘Se,f*nd  *5pect:  ratios  adopted 
devise  a  unified  theor?  whleh  f!ces-  11  is  desirable  to 

herent  in  nearly  all  of  *he  Inli  limit*tions  in- 

A  promising  avenue  in  this  direct  ion^entiJl?8  mentl?ned  above, 
character,  was  first  suBsea^Hk!  o?’  entirely  numerical  in 

21.  22)»  for  any  plinfoS  i{S  3ns^!"d.COli*bor,*tor3 
extended  to  planforms  «<«,  .2  al7  suP*rsonic  edges  (and  later 

arbitrary  def^t^n  «5e.  £hrtJC.l££!).0,J1“«1'«  *”  « 

of  the  "aerodynamic  influence  "8jb28ed  on  t*le  use 

pressure  developed  at  a  noiS?  In t£C5?nt'w  defined  as  «* 
velocity  of  the  fluid  (dSwnwash)  S  -If1!?  by  constant  normal 
wing,  while  the  downwash  is  IqlnmlU  elerentary  area  of  the 
planform.  This  influence  1coeffilSf»e*ro..i!VISr  the  reat  of  the 
adapted  to  loading  computation  ^iJ?te?ethod  can  readily  be 

laying  a  grid  of  boxes  on  £he  winged Slit  1m" Jh  by  °Ver' 
on  any  box  to  the  known  motions  of  In  Hi  J1  8  the  Pre8®ure 
applications  yield  quite  satisfaltoli1  th  boxe?-  Numerical 
forms  considered  in  Refs  21  ??  TM*CCUr!!Cj  ^or  bbe  Pi*n“ 

for  planforms  with  subsonic  ^dgei  b^ilt^  °?  wasuextended 
concept  of  Eward  (Ref  2^W  .£h  iy  introducing  the  diaphragm 

turbed  flow  region(off-the3)pUn?o?i*Cin8  b°XeS  °VCr  the  §ia' 

squares,  and^hi^can^e^hOSi0^6**!  tbe8e  references  are 
wings  with  subsonic  edges  to^thI°»-i Ehe  aPPlic*bility  for 
tabulations  for  the  aerodynamic  infi*®8  °£  ^  ^  •  Extensive 

sociated  with  souare  hollri^1  influence  coefficients  as- 

ing  •  suggestion  by  P^nes  LURefs*^1*^)  (Ref;  ?5)  •  follow- 
lar  analysis  using  as  ela^nt.l!,.3 ’  26  V  27 )  carried  out  simi- 
wbose  diagonals  are  peril Ul  to^h^Mlr^T/0"0*11**1  Mach  boxe8' 
M.cb  box  sppesrs  to  L^Lll 


•Reference  24  summarixes  the  essential  features 


of  Refs.  21-23. 
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nation,  Slojtld  by*L*\J!™P^utlng*inhu£!£e0f  ' f f ,  *PP rox  1 " 

another  alternative  grid  system  a*thors  suggested 

teristic  boxes.  These  el Up  of  s°“called  charac- 
sides  parallel  to  the  Mach^iSes*  ***  *r®  rho®bu,a*  with 

and  MortecchkyU(Refs?129^30),e for^nered”8?68****  by  Brandst«tter 
sonic  leading  and  trailing  »dff»a  pe*ed  planforras  with  super- 
harmonic  motions  Sey  ^SploS  ^bitrary  simple 

infinite  wing,  replacing^  J°iUt}°?  for  the  Wrter- 

one  ( the  foremost  points  of  thePtun  v**  a  rect*n8ular 

setting  the  downwash  equal  to  lewT^JE  h*1?8  the  8an,e)  *nd 
tangular  wing  between  the  the_re8ion  of  the  rec- 

This  step  is  pSSSJSbl •  si2j  3L*f2J40f  tWO  Planforms . 

The  downwash  is  represented  as  «h!Lif!din8<*d8?  13  supersonic, 
ates,  and  if  such  m  ?eS?esen?It?«£  ?! V®^6*  in  the  Gordin- 
must  take  a  large  numbEE  t0  be  8*tisfactory,  one 

tinuous  downwash  at  the  leadine  #di»e  *fC?knt  ^or  tbe  discon- 
The  method  of  Ref.  29  will  undfuht-fdi °fKtb®  taPered  planform. 
types  of  wing  performing  certain^osci?i EEt  useful  *>r  particular 
organized  in  such  a  wav  ^  lllations’  but  lt  is 

box  methods  from  the™£nputatioEEl  Efe*d  a?  efficient  a?  the 
types  of  region  must  bcTTreated  ™  standpoint.  Seven  distinct 
considered.  treated  on  the  most  general  planform 

that  they  ar^highly^EEtematized  thJ8®  num|rical  methods  is 

high-speed  machine c««pu2tlSl  d  ?he  letld e°re.weli  sui^d  to 
be  reduced  to  a  large  ageree**-®  «*TE®  affodynan,ic  theory  can 
erations,  which  involve*the8usp  repetitive,  elementary  op- 
tables.  However , ”0? epar ina  such^  f-K?ive58a11?  -gPlicAi; 
likely  to  be  an  expenEiEe  b3®S  for  the  flrst  time  is 

ficient  planning.  P  P  Ject  and  requires  careful  and  ef- 

recommend!JtioMtfor  inflienc^cSef ficient^hi 8  t0  Jurnish 

sent  a  set  of  working  rules  f«r°Ji!fiCient  tablea  and  to  pre¬ 
method  in  practical  flutter  analysesPr°The  h!®<°f  th®  box 
study  the  proposed  grid  system*  and  Tbe  baflc  problem  is  to 
various  tyfes  of  wln«  oFo^Lf  2  ?h!lr  *PP5  lotions  to 
sociated  iith  each  of  theso  »^?Lrl£uture  A.- 

culties  which  must  be  clearlvyr»^f„?h'5e  *re  diffi- 

procedureshfor'particular ^cases'6  “T  S5S3“' 

of  subsonic  edges,  must  also  be*tS2  S^KdiStS^KJ 
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«nd  investigated  seper«lli  fTS  !  probl“  have  been  isolated 
been  Studied  by  camarinl  th.  of  the  nethod  hal 

cording  to  the VJIi"  nS^ScIf  ^„f°r  *peclflc  ac- 

*  "  0th"  oerodynaaic theories ^f^knowr^valldityd* 

“°"*T»*dPin“tr«i"  cd.*»SLrtb:f*U  fjf  d't*11'd  investiga- 

2  jarsnas:  srs rss-T  srais- 

ssar Aasr 

tion  of  eerodJSini?  SflUe^c^  £oeffi°f  rul*8  for  th*  *PPlica- 
the  question  of  whet  dear^H  fficien£s»  one  is  faced  with 

ens wer  ^when°the  fp rob  lem^befnf^ena  &  to 

1  zpv, \?z  tUT.::  s  h‘:rn 

meted^lutte^speed^whereea^t*©^ 

very  insensitiv£.  ituELSaZ?**  «*•  speed  mJvb* 

sources  of  error  in  any  flSKl!*?!, to,  "^i**  «•«&  of  3! e 
Fhe  would  be  to  have  th*  *CU*at*on»  *nd  to  this  end 

loads  (compared  to  their  exarr  P?rcent*8e  error  of  the  air- 

Imown***’  k"*  order  of  'Mgnitude'SfiSrrc^T  values)  smaller  by 

srar-ss  ss  s.tssL? 

airloads  are  ;X?X*£in?in«°f  th* 

(1>  S!gh°!f £?i  *"<  so  that 

»ts.^^  ssJ.?32a«.*si“sr1. 

(2>  up  to  a  point 

ted  with  excessive  additional  ls  •«*°eU- 

in  accuracy^should  bt  ^ought^snd  “rked  Improvements 
when  they  yield  areatal and recommended,  but  onlv 
no5  interfere  seriously  with^h*10"*1  efficiency  and  do 
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SECTION  II 

AERODYNAMIC  INFLUENCE  COEFFICIENT  EXPRESSIONS 
FOR  THE  THREE  BASIC  TYPES  OF  GRID  SYSTEM 
II • 1  aerodynamic  Theory 

line?rized  Potential  flow  theory  for  simple 
harmonic  motion  of  an  almost-plane  lifting  surface  (Ref.  1?) 
the  complex  amplitude  of  the  velocity  potential  at  a 

point  fa#)  on  the  upper  side  of  the  surface  is*  ~ 

^Cx- 1)  _ 


where  arfap) 


Eq.  (2.1) 


is  the  complex  amplitude  of  the  normal 
component  of  fluid  velocity  (downwash) 
produced  by  the  wing's  motion, 


and 


Eq-  (2.2) 


are  rectangular  coordinates  (cf.  Fig.  II. 1), 

I®  the  deflection  amplitude  of  the  mean  sur¬ 
face  of  the  wing  at  (  %,  ?  )  (positive  down), 

A  is  the  entire  disturbed  region  bounded  by  the 

forward  Mach  lines  emanating  from  (  ), 

and  the  various  constants  are  defined  in  the  list  of  symbols/ 


Il!L?CCurdanc?  with  the  standard  procedure  for  representing 
simple  harmonic  motion,  the  actual  numerical  values  of  the 
potential  • ,  d^wash ,  etc . ,  are  the  real  parts  of 

the  SSiEr'SEHS?  etC"  Wh6re  "  18 
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*« iJpperhar,d°Tower  S&ef  ?«' 


f  *  h  ~k  =  zfb Lu)  +  Pjb Jffay) 


■  U)M1^  r} 

rAUfmw  u§‘  as(^-if-f^-f) 

deforaationSshapXeinastMdySf lowChth8  5l8S?ralc 'Po^omlal 

be  evaluated  i„Pt,r,  of  ?abul«ed  fS„J?Uble  cannot 

that  a  numerical  method  d  func^ons,  with  the  result 

idea  (Refs.  21-24??  one  Si  soUt^e^  t0;  /^^wing  Pines' 
small  elementary  areas  *-o  Illow  eertlil”*  A' i?to  sufficiently 
evaluation  of  the  double  integral  approximations  in  the 

downwash  is  constant  over  »ach  area  assVmes  that  the 

the  approximate  expression  ’  for  examPle>  one  obtains 

o 


HjtltT  tn^the  square'bracket|eof°EqWaf  g  17\  “!*  *  ■ 
influence  coefficient  (abbreviated  AIC  fr™  h*8  the  aerodynamic 
pressure  at  point  (  jcyJ  )  due  to  a  unf/S  herf  on)’  i*e*»  the 
The  AIC  will  depend  otf  the  W.n  i*  fownwa®h  over  A.  . 
and  the  shape  of  A  Three  kl!?n  °u  ^  relative  to  (  x,  v  ) 
for  the  elementary'areas  ( boxes?® ^ 10  shapes  have  been  proposed ' 
with  diagonals  parallel  to  the  Mark  ?^are  ?3xes’  ^tangles 

ss1  lelt°  Sacb MMs&si 
shapes  the  e  presJ^rio^rAjJIrfnoTd^^v'S888 
1 1  •  2  **““  J  * 


a. 
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EqU(2a4)d  which^re Snow°described  approxiinate  integrations  of 
Equation  (2.4)  may  be  put  into  the  form 

2SV?  ft#'  ] 


=  2sV2if(R+iZ) 


/ 


where 


£q.  (2.5) 


x>'a>  i-f‘  ■  z-f 

are  the  dimensionless  coordinates  /  rho  .  .  r  . 
square  box,  ’  ^  ls  the  side  of  the 

A,  =  "£  £ . 

'  Up* 


»/o'\  .  l^hL“~Ce«ndenv«fonsf  ^t'gr;tlon  correspond- 
■  »  Che  ouant  it  if.' »V*^°!?S  *_  **"  subscript  /  will 


be  dropped  froi" thTJSSKiS' ^TT  ’  T,  3^SCr‘pt  '  -ill 
ii'l,  the  pressure  at  the  center  Vk'  ’%*%:  As  shown  in  Fig. 
due  to  e  constenc  of-  -the  receiving  box  (n.m  ) 

dependent  on  the  distencesnWaSb  “  Che  sendln*  box  (  V,/J^  )  is 

Xn  x  -y  n  —  7/  r=  -2/ 

Wh"l  'box'  doTOStSe^^nd^be"^*"8;  ^  b£*  (n'^  >  ^  the 
arbitrarily  chosen  origin  box  f  cf°  theuright  of  an 

and  M  appear  within  the  brackets  of  E?  f^%?ar?roeterf 
v  and  ic  ,  the  eernriun.s'.T »  7  ci  or  ^  •  l2o)  along  with 

by  four  Independent  quIntmL  ^i!,C0f ff1?le?t  ls  £lxed 
coefficients  therefo?e\^e%  ■•faS3i£!.JSS5tJSi2.th— 
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Fig.  II. 1  The  Square  Box  Grid  System 

For  "far  boxes"  (i.e.,  v  >  Z  ),  the  following  types  of 
approximations  are  used: 

&  *'Cx-£c), 

'jf-#)  =  ^  67- %)=  C0S^' yfctf-tfrif) ,  Eqs  .  (  2 . 7a  -b ) 

where  (&,%,)  denotes  the  center  of  the  sending  box.  For 
sending  boxes  which  lie  completely  inside  the  forward  Mach 
lines  from  (  x,  ^  ),  the  areas  of  integration  are  squares. 

Then  one  obtains  for  the  AIC  the  expression 


R.J-il- 


mjp[-  (£  fakf-py)Kj 

rcos(v-i)l .  cos  (^vfp-ytpy) 

sm  ^4  •  ces(j(-/-F‘ -yy)  &**] 

+  .  os  (iyfe+£f-py)A%fi 

-  sm(v-i)l .  cos{ L 

-%,jjLCosv£,.cos(Z  B--J 
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Eq.  (2.8) 
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where 


v+-£ 


COS 


P(m+1) 
v  +  £ 


9 


v- 4 


cos 


-/  fi(/L  -t  /j  A  + 

V- X  * 


-  (*-i)Ai,p.  +efi-£)[‘Bt~'±L  -ati-'fdt  7 

ty-t)  efc-iy 


Eqs.  (2.9a-c) 


~  cc*h~' Eti  7 

@0*  +z)  @0*-  +z)  J 

fl5»tQ^ifndi?^b°X  t  ^ach  line»  the  «bove  expressions 

are  still  valid  provided  the  following  interpretations  are 
adopted: 


CiS~'x 

=  O 

fir 

Co&'x 

-  0 

% 

for 

0±X^/ 

L0S~'X 

“  JT~ 

f>r 

-1-SX.^O 

Cos'1*. 

-  7[ 

Xr 

x^-/ 

XX 

=  0 

for 

x±o 

cos/i'x 

-  0 

oc 

for 

COsh~'x 

=  0 

for 

-1  sx  S/ 

cosh  \ 

=  cas/T 

■'(-*) 

for 

-/ 

Eqs  . 

It  should  be  noted  that  only  positive  values  of  iZ 
considered  since  by  symmetry  ' 


( 2 . lOa-h) 
need  be 


/?_  _  +c  f _  _  —  ~  1~-c  I — 

A/ -A  Eq.  (2.11) 
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Equation  (2.8)  is  identical  with  Eqs.  (3a-b)  of  Ref  cs* 

?«Hnr  1  *  numerical  method  iffol-' 

lowed  in  Ref.  25  to  insure  sufficient  accuracy  at  hieh 

Sfvi£  by-.  &  ~?dtr*Jed  £,F1i;  «■*'  neat^boxissub- 

ence  length  Is  ^  A-  and  the  «ducedr£fequenc“i?'  ref'r" 


Sf  <M  O 


^i-8-  II. 2  Subdivision  of  Near  Boxes 

associated  with  the  small  boxes,  there  is  a  new  set  of  indices 
•  The  contributions  due  to  each  box  in  rows  *r'  =£ 
>T  /  “t7  thej  bl  c«lcul«ted  using  Eq.  (2.8),  employing 

*'  “  f//57  of  .  The  AIC  Sf  a  large  box^n  JoS8 

i?  obtained  by  summing  the  contributions  due  to  the 
corresponding  small  boxes.  For  example,  one  has  for  box 

2/=/J  M'-O  j 


7  ^ 


Eq.  (2.12) 


of  the'smaU^oxes  f^rows"  “”=7  «"Pnte£tbe^contr^utlo„s 

esC1“te  of  the  Practical  range  of  applicability 
of  this  theory,  these  rows  are  confined  to  the  large  box  ( a,  o  ) 
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!ni5?Ky^f0r"  *  triangular  region  A  as  shown  in  Fig  II  2  The 
contribution  of  the  whole  region  A  is  given  by  the  first  fJh 
terms  of  the  exact  infinite  series : 


Z7/A,f  sfo+izM *) 


Rj-iIA  =  )t 

A  A  ft  L  too  (  M  ' 


JZ8j  006 


<$>■] 


±r±/ A)_  9(**h) 

1? L/c  (^/  ljj)  J 


Eq-  (2.13) 


fro.  SS  corresponding  te^o^^f *£«•»« 

o”p“.!y  ^SSS*S.?JSi.uSg!»So?,h  not 
fi?ntheChighestnC/l  U  iS/f  ra°8t  fiv®  units  in  the  thi?drplaceS 

foT^cb-.Tber^r  (ie8r;s 

much  poorer,  as  will  be  shown  in  Section  III  4  q  (2'8'  1 
b.  Exact  Formulas 

funcM^?8  wf11"known  relations  in  the  theory  of  Bessel 

ver  ’Kyafkins  (private  communication;  no  formal  reference 

yet  available)  was  able  to  reduce  the  double  integral  of  £hp 
potential  function  [Eq.  (2.1)]  into  a  single  iSteSal  Sfth  an 
Jhfs8?^  exPr^s«d  in  terms  of  an  infinite  serief.  Following 

Sil  of Ea  3(?  lV  r?dUCtl°?  ls  P°ssible  for  the  double  inte- 
tabulated^ln  Rpf  ’inasmuch  as  the  AIC's  are  extensively 

II  2a ^th/pJ^'p25  USl?8  thf  approximate  formulas  of  Section 
«^2f*  the  exac5  expressions  for  the  square  grid  system  will 

E!  f.mPrS?ented  5eTe-  The  S«TC  "“thematical  technique  will 
slctlOTyiI13CUSSed  ln  connectlon  ulth  the  Mach  grid  systemln 
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where 


t  =  5?  '  Z  ‘  t 


.n.  ft 
/'  <$• 


*  *7 

are  the  dimensionless  coordinates  anx  4-  ^  ^ 

mens  ion  of  a  box  (the  spanvise  dimension  is*//*  S*reMtwise  dl‘ 

I  _  co4, 

'  vflL— 

r  = 

^J,i.  defineddd?f ferently3frorafl/^5  V‘“,hl.bost '  Here. 

0  ).  This  new  definition  ^  (by  a  factor  of 

flow,  the  AIC  becomes  independent  £or  3te"d!' 

a-  Approximate  Formulas 

similar  to^ha^fo^th^fqua?!  bSx  iEdT^d***  ejPressions  is 
here.  It  is  worth  notin^thlt  2  !?S  n??d,notube  repeated 
of  Mach  boxes  are  so  proportioned  thaf’i-Hi’?'  thed  tensions 
from  ( x.t  )  always  cut  the  boxes  aloia  ??rward  Mach  lines 
cuently,  for  complete  boxes  £ > Sa ^diagonal •  Conse- 
Chalf  boxes)  on  the  Mach  line  “?  f?r  Partial  ‘ 

symmetry  ^  /a/  •  Again,  because 

The  approximate  influence  coefficient  formulas  read: 

<**(*+ i)Ir  c*s(I  + 

i-cos(v-i)4r  cos(A 

-  j 


boxes 
of  the 


_  SL  _ 

-  i,p^.cos(g^t}£^  J 


(v*/zz>2) 

Eq-  (2.15) 
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(y,7})or(mt/i) 


where 


F18-  II -3  The  Mach  Box  Grid  System 


=  cas~'^z£  -  cos’' SlA 
v+z  v+£ 


y}~L  -  CCS~'  Si  i 

(2L  =  (v 


*~z 


-cos-'JOl  £ 


V- 


& 


JX.  (v  +*.) (*  £) yhp  -  -h(/i-~i)lasl~' _  cgsh V+A  7 

L  A~i 

■  (f*  +i)[ccsb  __  cosh~'  ?+  j  1  ,  -  _ 


Eqs.  f  2 . 16a-c) 
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quant itle*CvZ',reCat:^--n8  Eq*'  (2'}°«*>  »PPly.  The 

q  antities  and  G~~  are  independent  of  Mach 

%  s  **•  ^ 

in  Fig1Sir3“to  y?etd0f  112  b*  -P1^-  “’■•'own 

Eqs.  (2,17a-b) 

F°5  5uC  only  ^°x  (4*0  in  row  v=a  no  subdivision  is  necessarv 
and  the  contribution  of  this  box  ( Region  A)  is  necessary 


(Ra  +if  =  L/+  i/Af-£ 'MQ/ifl 


±f±/&, )_  (wiO/Af)  /-  -_„v 

M  l  z  I  AT/  6*77/  (yU^V-O) 


(^M 1-/) 

W  f 

Eq.  (2.18) 

Equation  (2.18)  is  obtained  from  the  infinite  series  expansion 
In  E$e  (“^a)XPreS8i°n  £°r  ^  «'-?•  •  whlch  *•  giver^below 

It  should  be  noted  here  that,  although  the  exact  expres- 
s ions  given  by  Li  for  the  AIC  are  correct  (cf.  Eq.  20,  page  18 
III'au’  £ertain  approximations  for  the  integrals,  n^cetfi-  ’ 
tated  by  the  type  of  computing  equipment  available,  are  in- 
adequate.  For  example,  Li  approximates  his  Eq.  (23b)  by  Eq 
(33).  Consider  the  exact  expression  for  the  steady-state  case 
with  constant  downwash  equal  to  unity: 

-  jJ  '{-(f-r-ity-iffo-tf-fr-if-]- ^ 

A--k  J 


*‘2 

V 


CoS 


tuL 

l-r+i 


Eq.  (2.19a) 

Li’s  approximate  formula  (Eq.  33)  for  4^s_f  ;  reads 
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The  integrand 
therefore  the 
singular  term 
firmed  by  the 
and  (2.19b) 


A  (f~A  (]-/*- iy'* 


Eq.  (2  19b) 

fta!?a  powlr1serufUi;Lf,the  up5er  «"<* 

leads  toScM*1*  ky  LI  for  the 

brief  table  beloS  This  fact  is  con- 

oeiow,  which  compares  Eqs.  (2.19a) 


reduce  the  dlmmslOTlIs^Ipnsslon'for^he’AlC,  ‘i.J?"11’1*  t0 

<*,),  v^RFFr~  ~Ul? 

Eq.  (2.20) 

ing^egioi!°boSSdednbvrtheClin«e(r  /"  A'  in  the  inf luenc- 

»slug3.  relation  whic&  can  be^rivVd  »„>• 

J"jo  (x  y IK  T>  )c*  a  TJT =  f,  x  ^ 

vxz-Ar 


v^Z2 
=  0 


for  X<& 

Eq.  (2.21) 
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Equation  (2.20)  may  be  put  in  thefonn 


±tiS.\  ±)  f* [\-3A-Q  CL  r _ , 

* '  ^'JJf  e  y#-( ifafEfJcxfoUzty 

,  Eq.  (2.22) 

Changing  the  order  of  inteeraMnn  anj  . 

P«rt,  one  has  gration,  and  carrying  out  the  7- 

i 

Eq.  (2.23) 

page  405  ofdRef°n32  relptions*  Eq-  ( 1) .  page  358  and  Eq.  (2), 

W*W*)  =  (a)Z(4)  +  2I](~0  Zr&Zrft) 

r^t 

roo 

-^smfcrsm'A)  fir 

Eqs.  ( 2 . 24a-b) 

and  since  !n/ .  for  thls  reglon>  ^ 

~t~/Z  ^~ZZr  (ffl  x~Z])jSl*(2rs/ri''#2£)  -snfcrsin  • 

^  7  TV  1  Eq.  (2.25) 

Finn wi-2  riuic1?}^  sv"8t*nt  md  the  in^ 

i:  ^  (&t  -  f  6£ftto  -  - j 


-  m.)-K&  . 

Equation  (2.25)  becomes 


Eq.  (2.26) 
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Eq.  (2.27) 


P.rtialI\o“StesV%durbiS2ch'f,f)1holds,JCrue  «"  f°r  the 
definitions  are  pbtmed:  Y  MaCh  Unes)  provided  the  following 


S/n  '%  -= 
S'n~'x=  ~ 
sir  'x  =  - 

z 


} 

J 


!  ^X 

-/  ±  K*J 
X  £  -/ 


Eqs.  (2.28a-c) 


and  A  unitsMawaybf^o®rJhe3poi«'f°r*r'‘)bOX5Wlth  ,lts  center  e 
The  following  expression's  " 


JL 

-2 


$££(«£)  (V-A-o) 


Eqs.  (2.29a-b) 
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where  £  is  the  familiar  function  of  linearized  unsteady 
supersonic  flow  theory  (cf.  Ref.  36  or  Section  III)  and  Jer 
is  the  Bessel  function  of  the  first  kind  of  order  Zr-  Al¬ 

though,  at  first  glance,  these  expressions  seem  rather  compli¬ 
cated,  the  evaluations  are  not  difficult  because  of  the  rapid 
convergence  of  the  infinite  series .  For  the  ranges  of  Mach 
number  and  reduced  frequency  of  interest,  only  the  first  three 
or  four  terms  need  be  retained.  Furthermore,  the  sine  functions 
can  be  reduced  to  Chebyshev  polynomials  Sn(x)  (Ref.  33)  when 
/x/  ^  / 


Stn  (zrsirf'x)  =(-/.’  A2  S£r,(  (2X) 
*  0 


1*1*1 

1*1  >/ 


Eq.  (2.30) 

where  S2r_.  (ZX)  are  represented  by  the  simple  expressions 


S,(ix)  =  zx 
S3(zx) -(zxf-eM 
-  fcxf- 

|  Eq.  (2.31) 
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c-  Steady- State  Formulas 

3teadi-n0wrp?obS'emsS\t1is  Mfch  b“  sche">«  for 

formulas  and  a  tabulation  for  the  AIC  "fot^I  corresPondlng 
(P.^a-b),  one  finds  all  the^l^  P.F«s  fq‘u.i  to  lTa  SB’ 

<&v  ’ -> 

*>.'  -  +T} 


Eqs.  ( 2 . 32a-d) 

The  values  of  are  given  in  Table  A  2  of  Ref  fnr- 


11,4  tn-yehafa^^IrlstloMrirg^^  M1“""  Coeff.o^e.  .a 
many  purposes  Sls1one'wlthtaxes°paranel*toCM»rhi?ite  8yste,n  for 

?ycJ  a  ""enslonless  system  wltE  origin  ac  ?“  *>"*«  ln°??eder 

trans format lonraCCerl3t^C  <  **  )  a/e  iefineWhe 


"ij&fc-Q-efri)] 


or  Eqs.  (2.33a-b) 

*-£-  &(*■**) 
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Eqs.  ( 2 . 33c-d) 


cnmnui 


Fig.  II. 4 


Characteristic  Coordinate  System 

“"‘ini  l*'*  M«e?oe.1”fttheo«st™tE3(Uif;1i’  tbe  ~lal 

gram  ( such  as  Region  A)  ii - —  fct  downwash  over  a  parallelo- 


where 


Eq.  (2.34) 


u)tfU 

Up 


Similarly  at  point  (  x -S,  #  ), 


the  potential  is 


Eq-  (2.35) 


Eq.  (2.36) 
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The  x  -derivative  of  if 
operation 


is  conveniently  obtained  from  the 


Eq.  (2.38) 

Similarly,  for  regions  B  and  C  of  Fig.  II. 5,  one  obtains 


Eq-  (239) 


Eq.  (2.40) 
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Fig-  II- 5  Basic  Regions  Under  Consideration  for  the 
Characteristic  System 


The  expressions  for  the  three  basic  regions  A  ,  8  ,  C 
may  be  combined  appropriately  to  yield  the  contribution  of  com¬ 
plete,  half  or  quarter -rhombic  boxes.  For  instance,  the  con¬ 
tribution  of  the  subregion  C (Fig.  II.  5)  is 


4,  fife  '5=0  -  k  fe  ~i '  i )  -fa  (rc -t.  sj  +f>  (rc-£  ,sc  -±) 


Eq.  (2.41) 


The  single  integrals  and  the  inner  portions  of  the  double 
integrals  in  Eqs.  ( 2 . 38) -( 2 .40)  can  be  expressed  explicitly  in 
terms  of  Fresnel  integrals  or  Bessel  functions.  For  instance, 
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From  the  definition  of  the  Fresnel 
tabulated  in  Ref.  34, 


Eq.  (2  42) 
integrals  which  have  been 


£{*)  =  C(x)-i  S(x.)  =  ~L  f 

Eq.  (2.43) 

one  may  deduce  the  following  identity 

^  #  ?(*v 

Eq.  (2.44) 

Therefore,  Eq.  (2.42)  becomes 


obtains81"8  reUtlons  sl”Uar  to  Eqs,  (2  21)  and  (2.2Ka-b).  one 


Eq.  (2.46) 


for2c Wo°hcYarious  fcyPes  obtainable 
(2. 45) -(2.46)  are  given  in  Ref.  28. 


Eqs.  ( 2 .38) -( 2.4o)  and 


a-  Steady-State  Formulas 


Since  frequent  use  will  be 
system  for  steady-flow  problems 


made  of  the  characteristic  grid 
,  the  corresponding  formulas  are 
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now  presented.  For  &  =  0  ,  Eqs.  (2.38-2.40)  reduce  to  the 

simple  expressions 


Eqs.  (2.47a-c) 

If  the  grid  numbering  system  is  that  of  Fig.  II. 6,  the  above 
expressions  yield,  respectively,  the  following  formulas  for 
the  complete  rhombus,  the  inverted  triangle  and  the  rieht-half 
rhombus : 


Fig.  Ii .6  Numbering  System  for  Characteristic  Boxes 
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i=  ££^/:a/ 

"  (3  -  'fl  L  X  1 


JU.  -  V  -  o 

^  =  0,v>0 


ZjUf  JL  /  _  *£/  1 1 

p  L  717  '"Jv+T  'JyjJ 

=  Afy/1  jJiggg  f  #+v  _  £2!lL-  ^tvH  I!  m,v>o 

f  L  ^  {'for, fan)  foui'  fo(t)ri)  \J^Lti)l>  JJ 


0) 


d) 


C"  -  r 

A,* 


f-  fC ■  fH 


■  fhi-jp*  ‘s-js 


Jj.=  v  =  o 

^0,  V>0 


■  5m 


p(2) 

AV 


c 


(Z) 

'A 


\-sm'IlAzL)]  «tv>o 

JU+V+/  /t+V-H  J J  f 


V^o 


cn\cm 

Av  VP 


cU)=  c 


Of 

P.v 


Eqs  .  (  2.48a -£  ) 

The  quantities  inside  the  square  brackets  of  Eqs.  (2.48a-f), 
»rePresent  the  aerodynamic  influence  coefficients. 
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SECTION  III 

APPLICATION  to  two-dimensional  problems 

coef  f  ic  lent  thSdi  "to  Cf  t  e  ywi!lJ  8  t  d?nan,i C  ' lnf  luenc  e  " 

study  various  aspect,  of  their  18  adv*nt««eou8  to 

flow  problems.  In  so  doing,  thePfill££f£2  t0  ^‘dimensional 
important  auestions  which  fin  be  nirifTi1?8  8re  three  of  the 
tive  simplicity:  *n  ba  partially  answered  with  rela- 

(1>  edge'b^.'Sroken  “  ”p"*?nt  ■  ««*Pt  leading 
is  n«ces.lt«ed  bv  mrf?.Ifpr0xl?!eion  BhlcS 
numerical  approach  ar.tJb^u^0"*  “ 

(2) 

<3)  chord*of  i°JJ*  distributed  along  the 

accuracy  is  afhievfd?  accePtable  degree  of 

can  also  study  thl^accuracy^f  the'tabflftT  the?e  purPoses»  one 
dynamic  influence  coefficients.  tabulations  for  the  aero- 


111,1  ^Presentation  of  Swept  Leadlm 


Broken  Linp« 


with  swetp*infleaAW°'dir:i<,ionai  8WePt  in  ateady  flow 

enough  so  that  the  parillel'TLdi™  th!i  Macb.nuniber  is  large 
supersonic.  If  the^ifg  pJodSffa  failing  edges  are8 

(constant  angle  of  attafkf,  tX  diSSIiSi  ' downwash  everywhere 
ference  at  each  point  is  constant;?!!? °nl®ss  Pressure  dif- 
formula  of  linearized  theory  (cf.  Ref.  given  by  the  exact 


(3  cot  A 
^l/Scot  5d  ~ , 


A  ^  ^  (3'^ 

”  /AS # ^the'square^nd*  the^Mach  «U  %L  *  *  ^  •  Since 
grid  system  la  placed  on  the  wlng®^,^*?’ 


-26-  • 

CONFIDENTIAL 


WADC  TR  56-97,  Part  1 


CONFIDENTIAL 


Mg.  in. 


S wep t UWing ^Us ingUMach 

/  -A  —  /)  1.  A  /  /  rrr-  I 


w  — — -  v  w*  u\|uai  c 

(&~0,  ta»A=  >/zj  M~  VE) 


£  sv¥  5s  °^*«:5My 

box«..by(M?kE  i^cZ?“y : 

^arsj^>tb?y2isi«ss ; 


h?vln* 

c5e  leading  edge  Is  approxi- 

thahhe  *id*®  °*  the  foremost 

Zf'lyZ'tJrS*  r**ort  to  this 
or  the  grid  system  is  such 

wing  taken  in  by  the  resultant 


cmonui 
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comma u 


ExperlenleUh«  ‘shorn  thJt°^hl?ewreaf  on,1“ed  the 
Ljl’-///?CUSCy  £°r  the  d«*™i"«i0n  of0il?loidt,1°nini<)'leld8 

reoeats  with  every  Second  col£X  \  f  h?jfWf8e  columns  of  boxes 
(3)  and  (5)  win  be  identicalTith  th^o  5  r,e?Vlts  for  columns 
fore  only  rows  (1)  and  (2)  need  hi  Il« <J°rJ 1]  ’  etc*  There- 
and  tabulations  of  the  AICs  (Table  Eq‘  (2-14) 

culate  the  approximate  dimension! pqq  nr-L  Re^‘  ^1),  one  can  ca.l- 
of  each  box.  The  results  for  b'  are^ho™*?  the  center 

cording  to  Eq.  (3.1),  the  exact  v*?»I  *5  ^  in  F1«*  IH.l.  Ac- 
points  is  '  exact  value  of  the  pressure  />'  at  all 


v  ^ 

and  the  dimensionless  lift  n@r  a*-* 

the  flow  is  per  unit  distance  perpendicular  to 


L/t  /un,/-  S/>  a/7 


(Sipi)Tf 


•  -=-/.  /sy  70 


Although  the  comparison  of  thl  individual6  Streani  direction, 
of  Fig.  m.i  wi?h  their  exact  vllue  SjI-?reMUre8  on  tha  »>oxes 
crepanci es,  the  average  lift  ne£  mile  18  aPPreci«Me  dis¬ 
close  to  the  exact.  P  Unit  apan  Proves  surprisingly 

the  wing  are  notrinfluenced6bySthePrr8°niC ’  the  Pressures  on 
ing  edge,  and  theref^fone  need  no!  fn!8iT*ti0n  °f  the  *«*!- 

sumed ‘cons tint  ^ver^ach 

suTiMr*®  sr-i.r^r A 

4*e*  o/  ffjZ / 

JJrtTi&t's  t^ntclon  of  the  p*rtui  »««•  fW/ 

(-/■  200,  5!  I  *)(**'*  off**/  ) 

4*e<z  ef  2 67/  / 

dealing /with#trailingaedgeai?ftPi6tIdb°r‘rf  this  wa?  of 

of  tb.  lift,  over  columns  (!)  Jnl'(Sf^S  &ffi/sh2Sr. 
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No.  of 

Chordwise 

Boxes 

J  Col.(l) 

ty) 

Col. (2) 

Average  Cols. 
(1)  and  (2) 

3 

'  -1.140,23 

-1-149.50 

-1.144,87 

4 

-1.143,83 

— 

-1.177,79 

-1.160,81 

5 

-1.161,75 

-1.133,58 

-1.147,67 

6 

-1.149,32 

-1-152,39 

-1.150,85 

Exact 

-1 . 154 ,70  | 

-1.154,70 

-1-154,70 

i«oie  ill, 


Exact  Lift  forL|[fTwoyDiron,erJCaliMet:hod8  with 
/  .  °r  *  Two-Dimensional  Swept  Wing, 

'  €-0,  tan -A.  =  '/z  ,  A/-  Vz ) 

fcuT  MS  lnc(ludes  «•••  «*«.  three 

the  chordwise  dimension  i?MUthCeeS^ly  *SSU“d  toefh[*Sp 

Severe!  conclusions  c.n  be  dr.wn  from  th„e  re.ult,; 

U)  becoraeVless°pronouncedSfor  boMr^Y*^ 

'&**'*> 

chordwise  ’boxesY^pracclcal0  ^ 
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(2)  For  calculating  the  lift,  which  is  the  iowe.t- 
order  generalized  force,  relatively  few  boxes 
However,  for  the  higher  order  gener- 
alized  forces,  such  as  the  " first  moment"  and 

ntL%Je?°nd  *h«  results  are  not  ex- 

chordwise  *°  '  *‘  *ccur*“  “lth 


from 


(3)  Since  the  pressures  fluctuate  appreciably  fi 

onl*  lo8ical  *nd  justifiable 
chordwise  integration  technique  is  the  rectangu- 

rule  assumes  that  the  pressure 
times  the  weighting  factor  associated  with  each 
J°de  J8  constant  over  each  box  and 
box*  *  t0  h*Ve  itS  value  at  the  center  of  the 

The  pressure  fluctuations  across  the  chord  ar« 
representation  of  the  leading  edge  The  tmasMoi  th 

«  <*1.  point  1.  whether  It  Is  possible  tt^iinp rove ^upo/th la* 

MrH«inK*tl00’  one  bad  th*  tabulations  for  the  AICs  for 
^oxes,  no  sucn  step  would  be  necessary  Hnu,uar  „  * 

tabulations  would  involve  not  only  the  swep  anglS  m  I  naJamll 
eter,  but  also  would  have  to  account  for  various  , 

of  the  boxes  by  the  leading  edgl  TheJcoJI  S 

atIlfirsrtoi>-^*t^0?  uf  prohibitive.  Another  alternative,  which 
c  “5ft18i*nce  n,^8ht  seem  reasonable,  is  to  take  the  AIC  for 

VP*"1*1  »  be  that  of  the  corresponding  c£pf«e  box 

times  an  area  factor  which  is  the  ratio  of  the  area  of  th* 
to  th*t  of  the  complete  box.  Upon  closer  stSdy 
tlin  proves  inadequate.  For  instance,  the  contribu- 

ls  c  ,th  Pressure/  at  point  <6  from  area  u7'c  of  Fig.  m  i 

Urgetr,; tZ7 'h* 

correction  had  becnused,  one  would  have  for 'the  c™tflb«”n 


(-  0.39Z)  (- 


af  a.  6  c 
A-zea  of  a6c 


)=-*■ 


/96 


as  compared  to  the  exact  value  of  -0.577  .  This  shows  that-  m 
simple  'area  correction"  for  the  determination  of  the  AIC  f^r 
partiai  boxes  is  not  satisfactory.  Incidentally  this  tvJ 
adjustment  is  suggested  by  Li  (pages  33-37,  Ref.  26)  ^Tf 
this  correction  had  been  used  for  the  determination  of  the 

results*8  **  d  ^  ^  in  the  flgUre  °ne  would^bt^n  thfpoorer 


£=-0.333,33. 
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t«ln  the  «*5e  of  «” |°ln|  .‘“"P1®  WlU  rc‘ 

£^r-^,3yS  *5  %^.£“ied 

t^e^fector’  X?™'  £S£  - 

•T/ £*,£/■"  ^^?S?r“lWsj& con- 

:» gSS  £  ?&.*&**  *nd1h*  ••~^?»3S2t:,  V 


Fig.  IH. 2 


SMDtUlnS1|lli!iT^a£l0?  ?"  •0Tw-Df*»»n*ion*l 
swept  wing  Using  Mach  (or  Square)  Boxes. 

f  ^  =  <?,  Af/J-4.  =  >/=  y^TJ 
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No.  of 

Chordwise 

Boxes 

Col.(l) 

Col. (2) 

Col. (3) 

. 

Average  Cols. 
(1)^||2)  and 

4 

-1.340,65 

-1.295,86 

-1.356,82 

-1.331,11 

5 

-1.329,04 

- - 

-1.335,60 

-1.331,58 

-1.332,07 

6 

-1.303,02 

-1.356,25 

-1.345,27 

-1.334,85 

Exact 

L 

-1.341,64 

-1.341,64 

-1.341,64 

t 

-1.341,64 

Table  III. 2 


?  Lift  by  Numerical  Methods  with 
Exact  Lift  for  a  Two-Dimensional  Swept  Wing. 

tost  A  =  Z/3  ,  m  -  Az) 


'.  =  0 . 


**  f 5* c  the  «ccuracies  of  the  lifts  for  individual 

columns  are  slightly  poorer  than  in  the  case  ficotA.  --  ^  but 

^;,*v*f*8*  Per  unit  span  over  the  three  foluwts  is  still 

being  within  1  per  cent.  Therefore,  when  1 

for  th*8rnrfi|0?^f  f°  Atf  ’  5°  Bi§ni*icant  change  in  accuracy 
for  the  total  lift  is  observed.  Of  course,  when  one  approaches 

(3  cot  A  -p- /  ,  the  accuracy  is  not  ex- 


- - —  w  a 

the  sonic  condition, 


.  ,  — '  .i — «tcur*cy  is  not  ex- 

pected  to  be  as  good,  but  is  probably  still  toleJable  for  wines 
of  finite  span.  (The  sonic  case  cannot  be  checked  with  a  two- 
w1?8'  The  above  and  other  examples  related  to 
two-dimensional  swept  wings  in  steady  flow  suaaeat  rh®  ....i,,,! 
th.t  th.  accuracy  1.  ....Stlally  “l?ora  fo?  mS  SlaS" 

?^Wtlon‘  yl?idln*  ?“tA  > '  ■  SiSci 

FlJ.  Ill  ?  .T??f  F  ?  ,y;t“‘  •"th»  •“».  the  result,  of 
fr8*/  I*1  „  J11-2  *lso  bear  on  the  case  of  the  square  arid 
A*  th*.1®^*r  Ma<:h  rofber,  M=/.z  ,  for  a  wing  with  sweep  8 
A  ~ -u„  (>/3)  ( (3 cd A  ^  /.  99) ,  the  results  for  the  square  arid 

are  shown  in  Fig.  III. 3  and  in  Table  III. 3.  Evidentl?  ?h?  ac- 
curacy  is  poorer  at  lower  Mach  number  when  square  boxes  are 
used  than  that  which  would  have  been  obtained  with  Mach  boxes. 
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Flg- 1113  *  T“°-DiMn*i<j"'1  ■-* 

(k  =  0,  tan  A  =  1/3  ,  M  =  /.&) 

$£• Sfi""?!  “lns  for  «htch  L  V5-’  “nd  I  *8  i"  v  /i“ 

b^nic^H  boundi”8  th®  boxes  are  arranged  so  that  the 

Fi«k  III1!  rinrrMltin8  th*  le?ding  *d|e  *•  »s  shown  in 
I  -'iZ11!4'  In  this  case  *  serious  difficulty  arises- 

Jalue  of  a*'*”””**™  *lwar  L°w  cohered  to  the  exact 

.  ,  ,  '  '  ^  .  Furthermore,  the  convereenca  t~ n — 

ward  thi3  value  is  very  slow  as  more  and  more  boxes  are  taken 
along  the  chord.  This  failure  may  be  explained  as  follows* 
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No.  of 

Chordwise 

Boxes 

ty) 

Col.(l) 

ty) 

Col. (2) 

ty) 

Col. (3) 

#») 

Average  Cols. 
(1)(2)  and  (3) 

4 

-1.087,26 

-1.170,43 

-1.095,73 

-1.117,81 

6 

-1.113,67 

-1.153,82 

-1.118,43 

-1.128,64 

Exact 

-1.156,65 

-1.156,65 

-1.156,65 

-1 . 156,65 

T«bl.  III. 3  Comp.ri.on  of  Lift  by  th.  Sou.ro  Grid  Sy.tom  with 
Exact  Lift  for  «  Two-Dimensional  Swept-Wing. 


Fig.  III. 4 


Pressure  Distribution  on  a  Two-Dimensional 
Wing  Using  Characteristic  Boxes. 

(> $=-0,  t<wA  =  '/Zj  m=  v<r ) 


Swept 
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effectively  repUc.d 

to P  "uTuSJ  MM?fv?“ely,?<,U‘',1'  th*  <=°n«ib<liion  of  «/c  * 

the  steedy-stete  ceef)  m*y,!’e  eeslly  proven  for 

tlv.  oueJity.  Siiuk  ^'"dw^oe'Io^p1^" 

*^y0"r  °"**^*'‘“°n’°f<>e‘ihi,chir«terl.[t!'g?S0'  '‘ 

•PP^enfly'repr^ent  thf*lefdin|  edge*byeiin^CpI^nn,l  °n5  **** 
pendicular  to  the  flow  as  in  t  lln®8  parallel  end  per- 

In  .o  doing,  on.  1.  f£«d  to  tibuUt£  the  ^C^i*^,?l8K^;-5)- 
boxe.  It  will  be  ehown  below  th.t  .«n  th.;.  hf??  t*lf‘rhomblc 

s^issa^sass. 


Fig.  III. 5 


Proper  Representation  of 
Characteristic  Grid  Systi 


the  Leading  Edge 
a  Is  Employed 


when 
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III. 2 


Adequacy  of  the  Assumption  of  Average  Con 
over  Each  Elementary  Area:  Steady'Hotion 


Constant  Downwash 


In  any  of  the  numerical  methods,  with  the  planforms  split 

d?f?ww?nt*ry  *r*?s'  j  basic  faaumption  is  that  the  downwash 
distribution  over  the  wing  can  be  adequately  approximated  by  a 

ff*  of  ?u}ta?ly  *ver*fe  constant  downwashes  over  these  areas. 
Although  it  is  possible  to  represent  the  downwash  in  each  ele- 
mentary  area  by  some  simple  function  other  than  a  constant  (such 
the  1  °5,  *h.e  aingular  downwash  associated  with  a  side 

edge,  which  will  be  discussed  below),  the  arbitrariness  of  these 
lunC5ionS  wo“ld  r«quire  an  undesirably  large  amount  of  tabulation. 
The  degree  of  precision  obtainable  using  the  constant -downwash 
assumption  is  closely  related  to  the  reduced  frequency  of  oscil- 
lation  and  to  the  rapidity  of  variation  of  local  angle  of  attack 
•long  the  chord,  as  will  now  be  shown. 

In  order  to  separate  the  effects  of  frequency  and  downwash 
variation  on  the  accuracy  of  the  results  for  pressure  distribu¬ 
tion  and  generalized  forces,  the  study  is  carried  out  in  two 
steps.  First,  a  straight  two-dimensional  wing  in  steady  flow  is 
considered  for  which  there  is  a  chordwise  variation  in  downwash. 
In  connection  with  this  point,  some  remarks  are  also  included 
for  the  case  when  there  is  variation  of  downwash  across  the  span. 
Secondly,  the  results  of  going  to  unsteady  flow  and  increasing 
the  reduced  frequency  are  examined. 

For  a  two-dimensional,  unswept  wing  in  steady  motion,  the 
pressure  is  a  point-function  of  tne  downwash,  i.e.. 


k- 


AT 


If) 


Cf 


Eq.  (3.2) 


The  pressure  calculated  by  AIC’s  at  the  center  l*.,# 

L  to 


Mach  box  J  n,m  rin'Fig. will  alao'*be>1’equal‘ To  the*  *"y 

quantity  {-■ur/o'  ),  because  all  rows  . .  contribute 

zero  to  the  pressure  y  at  (  *, ^) .  This  can  be  proved  easily, 
All  spanwlse  boxes  in  any  particular  row  ^  (  ir  /)  have 

the  same  downwash;  hence,  according  to  Eq.  (2.32d),  the  total 
contribution  of  this  row  to  the  pressure  at  )  is 


=  JL 

77 


Eq. 


(3.3) 
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Fig.  III. 6  A  Two-Dimensional,  Unswept  Wing  in  Steady  Flow. 
(Mach  Box  System) 


which  is  identically  zero  in  view  of  Eqs.  (2.17a-b). 
sure  at  (x,y  )  is  therefore  1 


The  pres- 


Eq-  (3.») 


A  similar  proof  holds  true  also  for  the  square  grid  system. 

Equation  (3-^)  implies  that,  for  the  two-dimensional  steady- 

state  case,  the  pressures  calculated  at  the  centers  of  Mach 

(or  square)  boxes  by  numerical  methods  are  exact,  regardless  of 

the  type  of  motion  in  the  chosdwise  direction  (provided  the 

constant  downwash  over  each  box  is  set  equal  to  the  value  of 

w  U*  ;enter> '  Specifically.  If  the  downwash 

x!  i?™  w1  eiln“f  l?  *  »  8ay  (  Ox/tt-  ,  the  pressure  *'  will 
be  given  by  the  relation  _  (^/^).  r 

If  two  chordwise  boxes  are  taken,  the  pressures  at  the 
centers  of  the  two  boxes  will  be  exactly 

fb=-o.2S  / 

where  ^  is  the  chord.  Using  the  rectangular  rule  (shown  in 
Sec.  II I. 1  to  be  the  numerical  integration  formula  which  one 
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should  adopt  in  general),  one  obtain,  for  the  lift 

I—  A#/*'?//  S/>*y 


4#)=.  W /“«’*■  SA 

9  TWJfw 


the  momentnaboutbthe ^leading  edge*'  How*v#r>  on«  obtain,  for 

fa  /,  )  _  M omen  f'/u*/  /  S/kxtn 


If,  say,  six 


✓  »>(#)**  -0.661,0  . 

tions  for  1th#'1|enerali£*d ' force*  *Ur*  *d*<*uat«  chordwise  integra- 
number  of  lnt.fr«"|  pllnS”?^.0"  *“*e  5*v*  *  •ufflcSl?” 
»ure.  .«  obt.ln.bl.  .wctly  «l?h^l!;%h!!?d2i.fl&“  th*  pr.,- 

Ion. tin t°f  1th*  c|“r»«erl.tlc 

’  the  P™»«ul.  MlSTunllll  *  ,lin“rly  "ith  x  , 
according  to  exact  theory  where* -  e»»I!**rly  *long  the  chord 

ch.r.cterls tic -grid  dcJutlirSIvSl.r*^”*^^*^^^  «» 
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iMi 


«chnIo;StonbePeJu.1?1to‘i1ttse™fu1r«C^  C?n8“nt  do'""“*h  °v« 

box,rC„e  has  for'the 


=  -  a.  zo¥,se 


/  '  =  -  0.38X34, 

'em  ' 


°i  r~/nd — • 

present  end  will  be  even  worss  f„In*  Y*fue,of  J*  is  always 
chord.  The  representation  of  Fig  ni^3*  *5£ker~*lon8  fc"e 
factory  whan  tfie  dowow.sh  v.rles^loiSg  ihe  cSrd  *  unS4Cls- 


box  equal  to  its^alue^at^he^nid  eho^'f^i!  d?wnwash  over  e«ch 
pressure  at  this  c^tr.J  polnt  FoJ  ^,fh',^1,*nd  ftnd  the 
down  washes  over  boxes  a  ..  ?  If  the  conRt*nt 

the  linearly  varying  case  as*’”  ^  'Fi8-  III. 8)  are  taken  from 


Fig.  III. 8 


®33»sM3Br- 


one  obtains  the  results  for  points  (/')  and  (  z'  ) 

fcf/>) ,oz> 2?  >  {i>)=  3o°j 

Tb...  contrast  with  the  corresponding  exect  value,  of 
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and  —  o.  3  7vS" 

tolerable. 


respectively.  Again,  the  error  Is  not 


Into  hal ^"rhombuses3 such  Jhat ° t ‘ f“U  rh°'”>“=  b°«* 
flow,  as  shown  In  Fig  m  f  ®  C~s  are  Perpendicular  to  the 


Flg-  111  9  Chirac  Cerlscl^Systernf  ln  ^  Edified 

are  boxes  ^  .  B, 

sures  thus  obtained  prove  to  bePexact  f„V  fr ’  »  the  Pres- 

Vfri5tions  when  the  motion  is  steady  if  typ?s  of  chordwise 

steady,  chordwise  strlnQ  /  Ql1^u  ^  the  inotion  is  un- 

make  contributions  to  the  pressS?e  a^t^oi^  J"  Fi8*  ni-9)  will 
be  shown  that  a  consistenter?o£  is  df  c’  Then*  lt  c*n 

*t  ,  since  the  downSashlllve?  Si  IS?*  f°r  the  Pressure 
those  at  the  downstream  edges  (such  as  are. taken  to  be 

average  values  at  the  centers  of  JvLff  (  *  rather  than 

trol  points  were  emploIS  ?  *  *  strips.  If  central  con- 

the  mid-chord  points  of  boxea  '  v"*  took  th«  downwash  at 

could  be  improved,  but  one  would"  then7  havi  11° 
for  an  odd™ shaped  area  such  as  tC  /I  Tf  . 

b,Jid  m  ....  ,,  ,uffuuntl,  ,MU 
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error  associated  with  taking  lower  points  rather  than  central 
points  as  control  stations  is  not  serious. 

When  a  three-dimensional  wing  has  swept  leading  edges  or 

parftllel*to  JhTJi"111  c?i8t  h*lf:rhon,blc  b°*«  with  on?  edge 
4  hf  Si?CfJthe  for«8°ing  analysis  has  shown 

«  desi5*bie  to  divide  the  rhombic  boxes  by  cuts  per¬ 

pendicular  to  the  flow,  quarter  rhombuses  will  anpiar  at  lulb 

boxes  *  ?h?«e<the  AIC  b*  tabulated  for  quarter  rhombic 

diaAd'wssntAg6,  since  the  number  of 

system*  1  the  t*bles  wil1  b®  four  tiroe8  th*t  for  a  Mach  grid 

tbi®  f*ct  in  view,  the  characteristic  grid  system  is 

"  S5l"b?2dm£S5  /PPlicttions.  Some  Imall  further 

use  will  be  made  of  the  system,  however,  in  special  cases  where 
no  chordwise  variation  of  downwash  exist. 

hnHpA‘  *  "®?c  instructive  example,  consider  linear  and  para- 

?Fii  TTTi?n^0n*pin  ?uwnw*fh  *®ro8S  the  span  of  a  finitewing 
IFig.  III.10).  For  the  sake  of  simplicity  at  this  stage. 


Fig.  III. io 


Notations  for  a  Finite  Rectangular  Wing  in 
Flow  with  Spanwise  Variation  in  Downwash 


Steady 


attention  is  focused  on  the  portion 
purely  supersonic  and  uninfluenced 
are  the  downwash  distributions 


A  of  the  planform 
by  the  wing  tip.  Of 


which  is 
interest 
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=  If  fay)  AT  fa)  =  Z/(/3yf 

unitlpan  £  rk^^TSJtai  f  “f«  P« 

actly  P  *  Keg  ion  f  fa#  -  ?  -  )  a*e  ex- 


a  anj  faft/«„,fsf«n  2  } 

(^0  *  (—*)  ~“u  r// 

downwash^the  dime™?™?^’  in  the  CaSe  of  llnear  variation  of 
,  e  dimfnsionless  pressure  is  constant  along  the 

ceSref  thP  ZqU*}*°-&  •  Usin8  any  the  nunwrSlpJo- 

n  'e  tr  !  T  pressures  and  the  lifts  per  unit  span 

?xacS  tthis  is  not  true  when  the  wing  has  a  swept  lead¬ 
ing  edge)  .  However,  when  the  downv;ash  distribution  is  para¬ 
bolic,  none  of  the  numerical  methods  yields  the  correctPresults 

lift  dLSSw?"  “I??8  fhe  'hord'  and  hence  the  2n“ia”5  “ 

aeaS  f^abfa°?r'?  a11  *1*°  d5Ylate  fr01”  the  <*««■  ",ay  be 
hv  -21?  J11*4’  these  discrepancies  can  be  alleviated 

by  taking  sufficient  number  of  boxes  along  the  chord. 


m 

Lift / Ur) if  Sfa.r> 

Lift /unit  S fan 

L'ft  / ur>,  f  Sfan 

t-f) 

L  Oiarfuj/se  MaJi  Lfoxes 

-  ■ 

(a  CLorfu/se.  fr/ach  fhres, 

(ff‘) 

fExacf 

0 

(1/6) 

( 1/^) 
(1/3) 
(1/2) 
(2/3) 
(3/4) 
(5/6) 

- - - 

-0.139,95 

-0.202,45 

-0.389,95 

-0.702,45 

-0.151,57 

-0.179,34 

-0.262,68 

-0.401,57 

-0.596,01 

-0.846,01 

-0 . 166 ,67 

-0.194,44 

-0.229,17 

-0.277,78 

-0.416,67 

-0.611,11 

-0.729,17 

-0.861,11 

Table  III. 4  Comparison  of  Lift  Distribution  According  to  the 

Mach-Box  Scheme  with  the  Exact  for  the  Purely  Super¬ 
sonic  Region  of  a  Rectangular  Wing  in  SteadyMotion 
(Downwash  Varying  Parabolically  Across  the  Span) 
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III. 3 


Adequacy  of  the  Assumption  of  Average  Const 
over  Eacn  Elementary  Area:  unsteady  Motion' 


Constant  Downwash 


v  In  Prevlous  section,  it  was  shown  that  the  square  and 
Mach  box  grid  systems  yield  the  exact  pressure  distribution  for 
a  two-dimensional  unswept  wing  in  steady  motion,  regardless  of 
the  variation  of  downwash  across  the  chord.  In  this  section 
the  same  problem  will  be  studied  in  unsteady  motion.  It  follows 

o£  fe!'  2  “5  E'l-  <J3)  chat  the  amplitude  of  pr«? 
sure  difference  between  the  upper  and  lower  surface  at  (  *  V  of 
a  two-dimensional  straight  oscillating  wing  is  1,0 


w  -  vit]  {-  *)  } 

%ihere  <ir  is  the  downwash  amplitude,  and 


Eq.  (3.5) 


I(u)= 


Eq.  (3.6) 


?fdth'f  ch^dh^Be?fCi  £uncelonuo£  <**  bind  of  order  zero, 

comes  Ch°rd  ia  t4ken  *s  the  reference  length,  Eq.  (3.5)  be- 


Eq.  (3.7) 


?ure  i  *re  dianens ionless  chordwise  coordinates,  A,  is 

the  reduced  frequency  based  on  the  semi-chord, 

is  the  modified  frequency  CD  =  zJkrf/p*  and  ^ 


Iff,*)  - 


Letting  ic -  %  =  u.  ,  one  has  for  Eq.  (3.7) 

-  E»-  (3-8) 


V 
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deflection  .?«y°po™t)°n  "U1  be  COnsldered  <  *  being  the 

(1)  2;  -  zSe‘Mt~  U,Jy  -  2i,k  ( vertical-translation 

oscillation) 

(2)  Z^Z-fae  (pitching  about  an  axis  at 

the  leading  edge) 


(3)  Zs-tife** 


_  Eqs.  (3.9a-c) 

•  li%=(2iAx*&]{  parabolic  chordwise  bending) 

For  the^irst  motion,  one  obtains  from  Eqs.  (3.8)  and  (3.9a), 

=  -2ti  {&C4x£  Cm,  Sx)-he‘MJt(^)\ 

7  Eq.  (3.10) 

where  f.(M,u ix)  is  the  xero-th  order  of  the  well-known  function 


i j.ii; 


£  as  •' 

Of  M  .  Similarly,  for  the  motions  (2)  and  (3),  one  has 

£=  x£(H hx )jn zx  - f-Kz^(Mtco<)l~  2k\-b  2tiJ 

+  x.%(M,ui*)\U*-}  ***•  (3.12a-b) 


*An  extensive  set  of  tabulations  for  the  A 
ly  being  prepared  by  the  Bureau  of  Standard 
wixl  appear  shortly  as  a  Technical  Note. 


-function  is  current - 
s  for  NACA,  and  it 
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fcWiiJS  bi*  dOVn”,h  ^  ““y  rather  than  2,4  , 


sq.  (3.13) 


guUrE££“na  io^and  Jo^ofVlS  *??f'f1alon.C?  .ch*  “»  lrl«- 

tribution  of  the  strip  A  l far  ^or  the  con- 

point  o  P  A  ( for  unit  )  to  the  pressure  at 


ti- e-*"d8fx,)J 

Eq.  (3.14) 


III. 11 


Illustrative  Example:  the  Use  of  the  MarK  . 
System  for  a  Two-Dimensional  Wi£a  wither  fee  4 d 
of  Downwash  Across  the  Chord  8  th  V*ri*tion 


If  th.  .trip  1.  taken  to  be  the  row  V  of  the  Mach  box  grid 
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mmam 


boxes  inE?hii3roi?  prSJldedtthl2||l?f  th*  Su®  of  tha  AIC  of  *11 

ITC  Mini  ava/1  •  **'  -  Cne _ ®X4CC  tXDriSfllftna  /  «  aa  % 


K-  Z  (#„+i£  )sy 

Jt=-V  *4  A  V,^‘L' 


Wh  K  Eq*  ^•15) 

chordwise  pressure  distribution*^  *°e  th?  ev*lu*tion  of  the 
one  gets  for  the  total  pressure  «%“«  *  "°tl°ns  (2)  and  <3>. 


so  Eq.  (?a6rT3r«rapp?JximtionWnWJshe*t  “id"?hord  of  strip  €7 
parameters  assume  the  values  "  exampl*>  let  the  ' 

&  =0-3,  0.99  W  ^ 

The  lowest  Mach  number  for  vMr>h  *.u  ,  , 

r.g.C..tef  C?  apply  satisfactorily  (  a/%  /%* theory  is  ex- 
omputations  since  conclusions  drawn  fnr  uS  ta^en  these 

conservatism.  ‘Wels^  of 

thio“h  ??iax”.CCOTpared  wlth  ««ct  resufts  L^gs'^nnr1' 

plots?eVeril  lnter«“"8  conclusions  say  be  dr.™  from  these 

boxes  ^yTeld'^satlsfactory^accuracj^for^the**  «  >'  a  f'“  chordwise 
cr«sV"  Also,  Pfor^motion™f!)Vthe*r  “f  "“”*«“»*  bo^s^'ln- 

than  those  for  motion  (2).  results  are  slightly  poorer 

b±8her  frequencies  (around  £  =  0  gQ\  ^ 

-ributions  using  six  boxes  aoa4n  ....  the  pressure  dis- 

number  of  boxes  wiU  yield  J5?J  -cceptable  but  a  lower 
not  shown  on  the  plot  for  tHe  LL  -  estiroa^ tea  (the  latter  are 
going  computation  was  madet  it  w^s  ffiSlSity^*  When  th«  fore- 
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- -  EXACT 

0  SIX  CHORDWISE  MACH  (or  SQUARE)  BOXES. 
EXACT  EXPRESSIONS  FOR  A  ICS 

x  THREE  CHORDWISE  MACH (or  SQUARE ) BOXES 
EXACT  EXPRESSIONS  FOR  AICS 


Fig.  III. 12  Real  Part  of  Chordwise  Pressure  Distribution  for  a 
Two-Dimensional  Unswept  Wing. 

(Motion:  Zi-^e^  &=0.3,  M  =  /.z) 
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-  EXACT 

o  SIX  CHORDWISE  MACH  (or  SQUARE)  BOXES 
EXACT  EXPRESSIONS  FOR  AlCs 

X  THREE  CHORDWISE  MACH  (or  SQUARE) BOXES 
EXACT  EXPRESSIONS  FOR  AlCs 


CHORDWISE  LOCATION 

2b 

Fig.  III. 13  Imaginary  Part  of  Chordwise  Pressure  Distribution  for 
a  Two-Dimensional  Unswept  Wing. 

(Motion:  Z=  zS-xe^  &=£>.*, 
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EXACT 


O  SIX  CHORDWISE  MACH  (or  SQUARE)  BOXES 
EXACT  EXPRESSIONS  FOR  AlCs 

X  THREE  CHORDWISE  MACK (or  SQUARE)  BOXES 
EXACT  EXPRESSIONS  FOR  AlCs 


CHORDWISE  LOCATION 

2b 

Fig.  III.  14.  Real  Part  of  Chordwise  Pressure  Distribution  for  a  Two 
Dimensional  Unswept  Wing  .  (Motion:  z  =  2Sxze^ fi=o.3  m ~/z) 
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EXACT 


O  SIX  CHORDWISE  MACH  (or  SQUARE)  BOXES. 
EXACT  EXPRESSIONS  FOR  AlCs 

X  THREE  CHORDWISE  MACH  (or  SQUARE)  BOXES 
EXACT  EXPRESSIONS  FOR  AlCs 


Fig.  III. 15.  Imaginary  Part  of  Chordwise  Pressure  Distribution 
for  a  Two-Dimensional  Unswept  Wing. 

(Motion:  Z~ 24xzeMty  -k^O.3, 
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REAL 

— A —  IMAGINARY 


EXACT 


* 

O  REAL  .  SIX  CHORDWISE  MACH  (or  SQUARE) 
X  IMAGINARY  BOXES.  EXACT  EXPRESSIONS  FOR 


CHORDWISE  LOCATION 


Fig.  III.  16  Chordwise  Pressure  Distribution  for  a  Two- 
Dimensional  Unswept  Wing.  (Motion: 
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-o-  REAL 
-A-  IMAGINARY 


EXACT 


O  REAL  .  SIX  CHORDWISE  MACH  (or  SQUARE)  BOXES 
X  IMAGINARY  EXACT  EXPRESSIONS  FOR  AIC*. 

4 


0  2  .4  .6  .8  1.0 


CHORDWISE  LOCATION 

ZD 

Fig.  III. 17  Chordwise  Pressure  Distribution  for  a  Two-Dimensional 
Unswept  Wing.  (Motion:  .*  =  £&&** 4*0.99  a 1  =  1.2.) 
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this  quantity  was  too  large  at  the  higher  value  of  k  .  Hence, 
calculations  were  carried  out  for  the  case  ty=/.s  ,  &=  <2.59  ; 
but  again  the  pressure  distributions  were  accurate  to  about  the 
same  degree  as  those  for  /*/*/•  2  ,  £  =  0.  99  .  It  has  there¬ 

fore  been  concluded  that  the  primary  parameter  is  it,  (which 
enters  the  formulas  through  the  expression  for  downwash)  rather 
than  .  As  &  increases,  it  appears  that  the  problem  requires 
more  and  more  chordwise  boxes,  regardless  of  the  Mach  number. 

A  more  meaningful  criterion  for  selecting  the  number  of 
boxes  is  furnished  by  a  comparison  of  the  generalized  forces 
calculated  by  any  of  the  numerical  schemes  with  their  exact 
values.  Using  the  rectangular  rule  for  the  chordwise  integra¬ 
tion  when  Mach  boxes  are  employed,  the  lift  and  pitching  moment 
about  the  leading  edge  per  unit  span  were  computed  for  three 
combinations  of  Mach  number  and  reduced  frequency.  These  are 
compared  with  the  exact  results  in  Tables  III.5-III.7. 


Z  Boxes 

3  Boxes 

G  Boxes 

Exact 

L 

k=0.3, 
M=1 . 2 

-0.204,96 

-0.352,52i 

-0.190,31 

-0.351,511 

-0.188,46 

-0.351,391 

k=0.99 

M=1 .2 

-0.378,43 

-0.645,69i 

-0.075,59 

-0.9l8,84i 

-0.064,91 

-0.923,601 

k=o.99 

M=1.5 

-0.070,43 

-1.221,551 

-0.046,62 

-1.247,611 

-0.044, 08 
-1.250,291 

k=0.3 
M=1 . 2 

-0.214,77 

-0.264,56i 

-0.202,87 

-0.250,811 

-0.200,73 

-0.249,321 

k-0.99 
M=1 . 2 

-0.259,38 

-0.741,941 

0.103,74 

-0.837,67i 

0.107,31 

-0.823,981 

(fy 

k=0.99 

M=1.5 

0.181,12 

-1.068,101 

0.171,80 
-1 . 106,431 

0.170,88 

-1.110,701 

Table  III. 5  Comparison  of  Lift  and  Moment  About  the  Leading 
Edge  per  Unit  Span,  Computed  Using  Mach  Boxes, 
with  the  Exact  Values  for  a  Two-Dimensional  Un¬ 
swept  Wing  Oscillating  in  the  Motion 
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2  Boxes 

3  Boxes 

6  Boxes 

Exac£~ 

4 

k=0. 3 
M=1 . 2 

-0.685,24 
+0.121, 36 i 

-0.678,70 

+0.102,551 

-0.674,72 

+0.091,08i 

-0.673,74 

+0.087,371 

k=0.99 
M=1 . 2 

-0.574,85 

-0.083,60i 

-0.591,53 

-0.461,831 

-0.583,42 

-0.478,831 

k=0 . 99 
M=1 .5 

- 1.029,47 
-0.l43,77i 

|  . 

-0.762,62 

-0.670,851 

-0.760,97 

-0.672,671 

n 

kO.3 
M=1 . 2 

-0.564,26 

+0.093,191 

-0.550,37 

+0.069,38i 

-0.542,43 

+0.054,851 

-0.540,07 

+0.047,931 

k=0 . 99 
M=l.  2 

-0.670,76 

-0.224,401 

-0.553,75 

-0.689,451 

-0.544,40 

-0.693,601 

k=0.99 
M=1 .5 

-1.085,55 

-0.169,421 

-0.711,45 

-0.962,88i 

-0.711,63 

-0.967,251 

Table  II 

1 - 

1.6  Comparison  of  Lift  and  Moment  About  the 

Hnit  Span’  ComPuted  Using  Mach 
with  the  Exact  Values  for  a  Two-Dimensi 
swept  Wing  Oscillating  in  the  Motion 
- - - 

Leading 
Boxes, 
onal  Un- 

-  zSxe iut 

2  Boxes 

3  Boxes 

6  Boxes 

Bxacjt 

L 

k=0.3 
M=1 . 2 

-0.793,41 

+0.171,561 

-0.801,59 
+0. 145,491 

-0.805,69 

+0.129,751 

-0.807,42 

+0.126,871 

Pf* 

k=0.99 
M=1 . 2 

-0.455,29 

+0.097,601 

-O.628.78 
-0 . 227, 83 i 

-0.622,44 

-0.264,051 

1 

k=0 . 99 
M=1.5 

-0.905,11 

-0.514,611 

-0.813,48 

-0.370,171 

-0.810,32 

-0.378,101 

M 

k=0. 3 
M=1 . 2 

-0.964,13 
+0 . 196, 40i 

-1.002,53 
+0 . 181 , OOi 

-1.024,60 
+0.171, Hi 

-1.036,06 
+0 . I67 ,97i 

a 

k=G.99 
M=1 . 2 

-0.606,44 

+0.082,721 

-0.814,23 

-0.386,47i 

-0.804,46 

-0.429,82i 

k=0.99 

M=1.5 

-1.146,76 

-0.796,561 

-1.049.32 

-0.580,761 

-1.050,65 

-0.593,841 

Table  III. 7 


Comparis°n  of  Lift  and  Moment  About  the  Leading 

Co,BPuted  Using  Mach  Boxes, 
with  the  Exact  Values  for  a  Two-Dimensional  Un¬ 
swept  Wing  Oscillating  in  the  Motion 
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The  following  points  msy  be  deduced: 

(1)  At  moderate  k's  (around  ^  =0.3)  and  for  all  three  motions 
where  the  pressure  distributions  are  satisfactory,  the  approxi¬ 
mate  lifts  and  moments  deviate  somewhat  from  the  exact,  but  im¬ 
prove  markedly  as  the  number  of  boxes  is  increased.  This  might 
be  expected  because  more  and  more  pressure  points  are  available 
for  the  chordwise  integration. 

(2)  At  higher  As  (around  A =0.99),  and  for  all  motions,  at 

least  six  chordwise  boxes  are  required  for  tolerable  accuracy. 
Tor  ,  the  integrated  lifts  and  moments  are  poorer  than 

those  for  M=/-5  at  the  same  A  This  fact  is  primarily 
due  to  an  insufficient  number  of  chordwise  points  for  integra¬ 
tion  of  the  sinuous  pressure  distributions  shown  in  Figs.  111.16 
and  III.  17.  For  A 1  = /.£  and  A  «•  0.99  ,  the  pressure  distri¬ 
butions  are  less  sinuous;  hence  one  obtains  the  better  accuracy. 

So  far  only  three  types  of  motion  and  the  two  lowest-order 
generalized  forces  have  been  considered.  If  higher-order 
generalized  forces  are  required,  which  is  often  the  case  in 
flutter  analysis,  more  than  six  chordwise  boxes  are  needed. 

The  same  requirement  must  be  imposed  if  motions  more  sinuous 
than  (2)  and  (3)  are  considered.  On  the  basis  of  assuming  that 
the  "half  wave"  of  each  such  chordwise  mode  is  adequately  repre¬ 
sented  by  a  second-order  polynomial,  one  may  extrapolate  and 
accept  the  number  6  as  a  minimum  for  chordwise  boxes  per  half 
wave  of  any  chordwise  mode. 

The  case  of  the  three-dimensional  wing  with  all  edges 
supersonic  (e.g.,  the  wide  delta)  represents  only  a  slight  ex¬ 
tension  of  the  case  of  the  two-dimensional  airfoil  discussed 
above.  Pines  et  al  (Refs.  21-24)  have  demonstrated  that  such 
wings  can  be  treated  with  satisfactory  accuracy  by  the  AIC 
methods.  Two  examples  of  this  type,  one  involving  steady 
motion  and  the  other  oscillations  at  a  moderate  reduced  fre¬ 
quency,  are  presented  in  Section  V.  From  examination  of  those 
calculations,  in  conjunction  with  the  vabious  two-dimensional 
cases  given  above,  it  has  been  concluded  that  a  conservative 
rule  for  wings  with  all  supersonic  edges  is  that  a  minimum  of 
eight  boxes  should  be  taken  along  the  midspan  chord.  This 
minimum  should  be  increased  if  the  wavelength  of  the  mode  of 
chordwise  deformation  becomes  too  short. 

The  rule  involves  a  slight  upward  revision  of  the  re¬ 
quirement  of  six  chordwise  boxes  recommended  for  purely  two- 
dimensional  cases.  This  modification  is  introduced  to  account 
for  the  possibility  of  added  errors  which  may  arise  due  to 
spanwise  deformation  of  the  finite  wing,  since  these  errors 
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may  be  additive  upon  those  due  to  chordwise  deformation  and 
already  reduced  to  acceptable  size  by  the  six -box  rule.  Another 
consideration  affecting  the  decision  to  recormnend  eight  boxes 
is  that  the  increment  in  computational  labor  caused  by  going 
from  six  to  eight  is  not  severe. 


IH-*  Comments  on  the  Accuracy  of  the  Approximate  Formulas  for 
the  Aerodynamic  Influence  Coefficients. 

The  approximate  formulas  given  in  Section  II  for  the 
square  and  Mach-box  grid  systems  nwst  be  investigated  for  their 
precision  at  high  values  of  &  .  The  most  extreme  case  in 

the  tabulations  of  Ref.  25  (4=**,  will  be  tested.  For 

rows  v  -  o  and  /  ,  some  caution  was  exercised  in  Ref.  25  for 

high  values  of  6  ;  a  subdivision  technique  which  splits  the 
basic  square  area  into  25  smaller  squares  was  introduced  to 
improve  the  accuracy.  However,  for  rows  v>2  ,  no  such  step 

was  taken,  because  it  would  entail  a  much  lengthier  set  of  cal¬ 
culations  .  One  can  determine  the  order  of  accuracy  for  tT~  z 
=  and  by  using  a  subdivision  technique.  For  in-’ 

5be  Urge  box  *=z  (Fi«-  H  I)  is  taken  and 

subdivided  into  twenty-five  smaller  boxes,  one  has 


Eq.  (3.17) 


Since  tabulations  at  this  Mach  number  are  also  available  from 
Ref.  25  for  A^cos  ,  this  check  can  be  carried  out  very  simply. 
Table  III. 8  shows  the  comparison  of  the  AIC's  for  , 

2  given  in  Ref.  25  with  those  calculated  by 
Eq.  (3-17).  J 


It  is  seen  that  all  entries  in  the  table  are  appreciably 
different  from  those  computed  by  the  subdivision  method.  The 
latter  set,  being  obtained  from  a  finer  grid,  must  necessaril' 
be  much  closer  to  the  exact  values  (this  will  be  proved  later' 
A  similar  calculation  could  not  be  carried  out  for  all  boxes 
(v^3)  ,  there  being  insufficient  entries  in  the  tables  of 
Ref.  25.  Although  a  direct  comparison  of  the  values  in  the 
first  column  of  Table  III. 8  with  the  corresponding  exact 
values  would  have  been  preferable,  such  an  undertaking  would 
have  required  the  difficult  task  of  evaluating  the  AIC's  from 
the  exact  formulas  for  a  square  box. 
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From  Ref.  25,^=44' 
M=/.Z 

With  Subdivision,  Using 
Ref .  25 ,  ^  =  a.  oe,  a i=/.  z 

%  +  L 

-0.091,77-0.004,241 

-0.099,30-0.003, 14 i 

Rzj  +  L  ^~2j> 

-0.107,66-0.036,251 

-0.117,82-0.037,881 

R2/Z-*-lI2Z 

-0.210,49-0. 262, 8li 

-0.236,11-0.258,861 

-0.060,88+0.077,151 

-0.020,95+0.045,891 

0.176,40-0.016,241 

0.153, 12-0. 011, 67i 

Table  III. 8  Aerodynamic  Influence  Coefficients  in  Row 
-v-Z  Calculated  with  and  without  Sub¬ 
division.  =  o.^M^t.2.) 


„  »?th  than  COBnP*re  the  accuracy  for  the  individual  boxes 
in  a  given  row,  one  can  more  easily  study  the  sum  of  the  AIC  in 
^urow-  As  was  pointed  out  in  Section  III. 3,  the  exact  sum 

(3  if)  C°TiblfUTTTnQ  b°XeS  ^  *  r°W  iS  glven  ^  E^* 

Table  III. 9  shows  this  sum,  computed  in  three  dif-  * 

ferent  ways,  and  indicates  significant  deviations  for  rows 

in  SiS  value  of  4'  *nd  M-  With  subdivisions 

in  row  -u—z  ,  the  accuracy  improves  markedly. 

Additional  calculations  for  lower  values  of  at  this 

shorn  Jor  the  approximate  fomulas,  the 

maximum  tolerable  value  for  is  around  o.oe  if  one  is  to 

obtain  adequate  estimations  of  ^  •  At  higher  Mach  numbers, 

0Ahm  f!34f?alncr^ased  somewhat ,  but  in  all  cases  the  value 
ot  the  modified  reduced  frequency  should  be  lower  than 
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^  0.3 


Eq.  (3.16) 


V 

Ref  .25,  =  0.  ¥J  M-/.Z 

With  Subdivisions 
(Table  III. 8) 

Exact  [Eq.  (3.14)] 

K 

0 

-0.794,33+0.382,331 

-0.796,26+0.382,621 

1 

+0.650,46-0.254,421 

+0.651,92-0.251,391 

2 

‘0.329,69-0.318,751 

-0.360,07-0.350,351 

-0.364,61-0.352,961 

3 

-0.021,63+0.312,311 

+0.009, 21+0. 302, 53i 

4 

0.116, 31-0. 307, Oli 

0.131,36-0.265,071 
- 1 

Table  III. 9 


The  Sum  of  the  Aerodynamic  Influence  Coefficients 
in  Rows  Z=c,  ...  ,  ^  for  4/  =  0.V, 


M=f»S.  ,  Calculated  by  Various  Methods 


It  should  be  pointed  out  that  the  majority  of  the  tables  in  Ref 

%:?!*£  ssi.sss's-ruigj.s  Mbu!.«oSct1!?  “„3b  be 
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tables  are  to  be  based  on  the  Mach-box  scheme,  no  serious  diffi 

high-speed ^m,chinePcomput«t ion. *^*PCln^  to 

of  Re^f0"^"^^ 

cp-Sf'ni^c 

tweJn  th«  LJ2SrS1£5K°f4the  receiyin8  point"  but  still  be- 

Llh".s  20  ciS?Lis2e^“™bfb^”eJ«S4p^  i„I”S' 

?hil  tac^ie™n«s°nvenlenc*  C*n  be  ““’“ted  to  result  free 
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SECTION  IV 

LOADING  OF  PLANFORMS  WITH  SUBSONIC  EDGES 

the  wide  delta%ing)PUthlynrpPerS°niC  surfaces  (e.g., 

influenced  by  any  area  0f f^h^nf63*00  P^an^orni  are  not 
planfonn  has' subsonic  edges  ther^f0?1^  .?owever»  when  the 
the  x-u-  -plane  adjacent  to  re«ions  of 

sure  distribution  on  the  wholp  nr  w*1*c:*1  affect  the  pres¬ 
lowing  Eward : s  concept  (Ref  S“ff S«-?5^ the  Planforn»-  Fol- 
permeable  'Miaohragm"™!, tl \  3?)  5  fi<;titious,  thin,  in.- 

diaphragm  is  aLumfd  to  coincide *2? th*!!  these, regions .  This 
therefore  it  cannot  al?e£rhf??  ith  a  Strearn  sheet»  and 

sure  difference™?^  P^i«  bet™ennusC?Llt  ??“ln  *  P«- 
faces.  The  corabin*tJnn  tween  its  top  and  bottom  sur- 

a  new  surface,  which  is  purely  ^ers^ic^nd^  forn,s 

washes  are  known  on  the  planform  P«h  for  v*ich  down- 

zero  over  the  remaining  area  tAp  Jjf1®  the  Press“re  jump  is 
strengths  of  sources  and  si^«,  proJle?  is  to  find  th£ 
distribution)  which  must  be  plaLd  the  downwash 

then  be  d.t^'2  m.y 

used,"?«U£«n™CJf  as  Che  AIC  method  i,  to  be 

diaphragm  region  is  necessarilvPlim??prf  JondiJJ°"  on  the 
points  of  this  region  A^  nn,‘  mf^ed  to  a  finite  number  of 

in  general  be  zerT  The  Ire wil1 
the  effects  of  this  appwxiMtio£  The  ??  for 

number  of  boxes  on  the  di-nhJ™  Jh  first  is  to  take  a  large 
boxes  but  «?suL°Eef??eh»§hr^w«b  ’trutlt,*  tC  taka  £e“« 
cWon^tiln3PeCUUy  adJ“«ed  *>  egree^ith^he ^zero-pressure**0* 


slbleILrtods“^rt?e;l«?p1«  ?rfe«lrC'l0nKt0  fug«ese  P»- 

flrst,  the  p«rtlcul?M£,?”S?«  °%  ,ubs?nlc  edges.  At 

!Si nlil.?5  ssjjsj  «-r'd°^solfce 

cases°of  £.dy 

™SSS£S°SJE5E: for  “iSS’SS  &J°E  Ss. 
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IV .  1  The  Side-Edge  Problem  In  Steady  Flow 

Consider  a  rectangular  wing  in  steady  flow  at  a  constant 
angle  of  attack  <x  (Fig.  IV. 1).  The  aspect  ratio  is  sufficient¬ 
ly  large  that  the  Mach  lines  at-  and  cd  intersect  behind  the 
trailing  edge  (this  is  not  a  necessary  restriction  on  the 
method  of  approach) .  The  exact  dimensionless  pressure  dif¬ 
ference  at  any  point  on  the  "mixed"  region  B.  is  given  by 
Lighthill  in  Ref.  39  to  be 


Eq.  (4.1) 


At  any  point  in  the  purely  supersonic  region  A  ,  one  has  the 
Ackeret  result 


A'a  At  =  _  1 

W)' 


Eq.  (4.2) 
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Because  of  the  manner  in  which  x3=  , 

Eqs.  (4.1)  and  (4.2)  dimensionless  £'  '  jLeinpl5yed  t0  fflake 

srsks/s  gi.'Sbissa?  c,icul*ted  by  ■ssfe'.sSi. 

ing  me thods°*the Character 1st ic°gride  l0id‘ 

5&" 

Section  III)  is  eliminated  It  al-sn^f  fUS6  thJs  systern  (see 

being  able  to  satisfv  t-ho  offers  the  advantage  of 

points  on  the  side  edge,  a  fac^which^llf"^1610"  at  dlscrete 
the  procedure  is  illustrated  *PP*rent  «s 

boxes  ere  assu*d,  ..  ^ 


Flg.  rr..  ^»-t5ort1^„1'.?rc“  sSnSr4°Ieao?e^J"f£lar  Wlns  ln 

sumptions : CW°  d:fferenC  w,),s’  based  on  the  following  two  as- 

<l>  Slv»;&s:  s’si&s: 

the  downwash)  is  assumed  constant. 
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°j  ss  e°3gethe  ssfnsjB1- 

issx  zgsgms  %~k  }nn3jrI  T  °*ution 

proportional! tv  i.  I*  *nd  involves  inverse 

distance  off  t^edge.  qUare  r°ot  of  Perpendicular 

Procedure  1; 

rrr*”  r  *■  ^  s-itfe-w*  - 

the  relation  *  condition  of  zero  pressure  at  point  0)  yields 


or 


Eq.  (4.3*) 


(C+  c:‘) 


where  the  AIC's  CA  are  those  of  Eqs.  (2.48) 
d  it  ion  for  p^int  (z.)  leads  to 


Eq.  (4.3b) 
The  same  con- 


co^-  - 


O) 


( S'  7 


A-5-1  ^  (z) 

C  ^  4a 


) 


Eq-  (4.4) 


strengths  "““‘V  S°luti?n  for  the  down  wash 

a  r6j  successively  ejuai  to  i2£ln8f£!!  pr®ssu^es  « 
need  not  be  considered  since  theS  5?ro*  '5°te  th«t  boxes  ?-io 

at  the  discrete  potats’of  fate«s?  )  "°C/“ect  “e  pressures 
•re  determined,  the  pressure  at  «  ?  fe  these  downwashes 

calculated  by  simple  summation  of  the°effects°of 'all^bSS  Ca?  be 
the  forward  Mach  cone  emanating  fromST  The  5a?U  bo£ei,ln 
are  shown  in  Table  IV  1  and  Sa.il  Pe  values  of  A 

to  Eq.  (4.1).  Thl  accuracy  3£S?J2!i^Jh  the  exact  according 
factory.  As  the  number  of^oxes  ifb}e  iS  n°5  re8arded  as  satis- 
prove  somewhat ,  huHE  «nv£«»«  S  '.!he  re|ult?  in¬ 
stance.  if  eight  boxes  had  be^en^^d  Se^u/V"’ 
would  have  the  value  of  -o.z73,t.i  (  ^  wtch  £our  bo^s) 


-63- 

CONHDENTIAL 


WADC  TR  56 -97 ,  Part  1 


C0NF1DENTUI 


AsS£ill  be  seen^ater  this^oSv^^  t0.the  exacC  ~a  3*/  83 
or  square  boxes  are  employed.  The^easo^for^h^S^  When  Mach 
xn  the  present  calculation  is  that-  f  5he  discrepancies 

. diaPhragm  box 


Procedure  (l) 

Procedure  (2) 
Downwash 
Singularity 
Included 

••  ~ requir 

Exact  [Eq.  (4.1)] 

| 

P12 

-0.254,62 

-0.401,07 

-0.391,83 

H, 

-0.171,97 

-0.291,31 

! - 1 - - 

-0.295,17 

P19 

-0.135,03 

-0.238,01 

-0.246,75 

p14 

-0.377,99 

-0.507,75 

-0.5 

Pis 

-0.273,62 

-0*389,33 

-0.391,83 

p23 

-0.221,55 

-0.326,29 

-0.333.33 

pi? 

-0.454,83 

-0.570,86 

-0.564,09 

p22 

-0.334,10 

-0.452,62 

-0.454,37 

P21 

-0.508,67 

-0.614,25 

-0.608,17 

p26 

-0.396,95 

-0.498,70 

-0.5 

>k 

-0.549,12 

-0.646,54 

-0.640,98 

p28 

-0.580,96 

-0.671,82 

-0.666,67 

Table  IV. 1 


:S!“S“nS«SSi:“' 

B™5«S'dU‘  ^  F°™’  of  lh"  Ch.r.«^«ic 
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,  Ifc  is  known  from  Lagerstroro' s  exact  solution  (Ref.  40)  that 
the  downwash  over  the  diaphragm  region  exhibits  a  square-root 

Jk  fc5e  Sid^  is  aPP5oached  from  theqoutside  and 
vanishes  on  the  forward  Mach  line  ^  .  The  effect  of  this 

singularity  can  be  approximately  accounted  for.  In  Cartesian 
coordinates,  the  receiving  point  is  while  ?  v  are  the 

IV™??"8  A^h^16?  r®pr®sentin8  the  influencing  point?(  see  Fig. 
i  iA  ^h?TaC5®ri3^iCu^ coordinate  system  with  origin  at 

*s  defined  such  that  the  axes  ^  and  s  are  resnectlv- 
ely,  the  right  and  left  forward  Mcch  lines  emanating’ from 


r-Jp Sfr-rJ] 

s- 


Eqs .  (4 . 5a-b) 


‘S'  r“elvlng  P°JRt  on  the  side  edge  where 

dtaphragra°region  £°r  do"""*sh  distribution  on  the 


4J~(rJs)  ~  zsl 


K-r 


Vd  Vr-s 


(  r>s) 


Eq.  (4.6) 


Equation  (4.6)  exhibits  the  proper  singular  behavior  at  the 
edge  and  vanishes  on  the  Mach  line  -t£-c  (  r=  ^  ).  Here  d  is  , 
typical  length**  and  is  the  "strength"  of  the  downwash.  If 

the  receiving  point  is  on  the  diaphragm  region,  Eq.  (4.6)  be¬ 
comes  ' 


-tirU,  s)=  ~  r  (Fy.iv.va.) 

and  if  the  receiving  point  is  on  the  planform, 


Eq.  (4.7a) 


t5is^!ectj2n  th®  coordinates  are  dimensional 
in  contrast  to  Section  II. 


**For  convenience,  this  typical 
derivations  f  u.  the  side  of  a 
a  Mach  box. 


length  cL  is  taken  in  subsequent 
rhombus  or  half  the  diagonal  of 
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Eq.  (4.7b) 


^7^5^=  ykrJlZl L__ 

'^(r-r0 )-s 


vrtiere  the  notation  is  that  of  Figs.  IV. 4. 


Fig.  IV. 3 
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One  can  easily  prove  by  direct  substitution  of  Eqs.  (4.7a-b)  into 
the  characteristic -coordinate  form  of  Eq.  (2.3)  that  the  as¬ 
sociated  pressure  influence  coefficients  are  given  by  (for  <o  =  o  ) 

_  fK  - 


c  - 


(tfs) 


fir  r,j 


/7*.  MV*. 


<f  fir  Fy.  /**'. 


Eqs.  (4.8a-b) 


Returning  to  Fig.  IV. 2,  the  following  downwash  distribu¬ 
tions  are  assumed: 


(1)  A  singular  downwash  distribution  over  the 
region  bounded  by  .  Additional 

singular  distributions  uj^i  ,  over  the 

regions  a-OU  and  a(3)e  ,  respectively. 

(2)  Additional  constant  downwashes  and 

over  boxes  Z  ,  V  and  S'  . 

One  obtains  from  the  zero -pres sure  condition  at  points  0)-(6) 


=o 


Of 

«£  Co,o 


Eq.  (4.9) 


as  in  Procedure  1  a  sequential  solution  yields  the  required 
values  for  the  assumed  downwashes.  Incidentally,  this  sequential 
solution  is  equivalent  to  the  inversion  of  a  triangular  matrix 
as  shown  by  Pines  and  Dugundji  in  Ref.  24,  if  the  equations  are 
cast  in  matrix  form.  Once  ,  are  determined,  the 

pressure  at  any  point  in  the  mixed  wing  region  can  be  calculated 
For  Instance, 
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Eq.  (4.10) 


The  pressures  obtained  in  this  manner  are  also  shown  In 
Table  IV. 1.  A  marked  Improvement  Is  evidently  achieved  when 
the  singularity  of  downwash  Is  accounted  for,  even  though  the 
zero-pressure  condition  Is  satisfied  at  the  same  number  of 
points  on  the  diaphragm  as  in  Procedure  (1). 


Next  considered  is  the  Mach  grid  system  for  the  same  wing, 
with  six  chordwise  boxes.  As  may  be  seen  from  Fig.  IV. 5,  two 
difficulties  arise: 

(1)  The  pressures  cannot  be  made  identically  zero  at  points  on 
the  side  edge  since  the  control  points  (centers  of  the  nearest 
diaphragm  boxes)  lie  off  the  wing; 

(2)  if  Procedure  (1)  above  is  to  be  followed,  no  disturbances 
can  be  assumed  in  the  triangular  regions  -Z>  unless  additional 
tables  for  such  regions  are  prepared. 


Fig.  IV. 5  Mach  Grid  System  for  a  Rectangular  Wing  in  Steady 
Motion  at  Constant  Angle  of  Attack 
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Table  IV. 2  presents  the  pressures  at  the  centers  of  t-h*. 
*  calculated  accordine  to  t  n _ j  /~\S.  t"e 


■ 

Procedure  (l) 

Procedure  (2) 
Downvash 
Singularity 
Included 

Exact 

t Eq.  (4.1)] 

pll 

-0.608,17 

-0.445,87 

-0.391,83 

P12 

-0.410,57 

-0.326,22 

-0.295,17 

p13 

-0.352,90 

-0.248,31 

-0.246,75 

pU 

-0.300,01 

-0.233,31 

-0.216,35 

pi5 

-0.274,05 

-0.191,69 

-0.194,98 

p18 

-0.704,83 

-0.599,84 

-0 .564,09 

P19 

-0.525,33 

-0.475,17 

-0.454,37 

P20 

-0.460,54 

-0.391,79 

-0.391,83 

P21 

-0.401,99 

-0.362,09 

-0.349,80 

P25 

-0.753,25 

-0.669,62 

-0.640,98 

p26 

-0.590,81 

-0.552,19 

-0.535,44 

P27 

-0.525,88 

-0.470,63 

-0.471,02 

P32 

-0.783,65 

-0.712,07 

-0.687,49 

p33 

-0.634,79 

-0.602,35 

-0.587,93 

p39 

-0.805,02 

-0.741,43 

-0.719,56 

Table  IV. 2  Dimensionless  Pressures  on  the  Mixed  Region  of  a 
Rectangular  Wing  in  Steady  Supersonic  Flow  at  a 
Constant  Angle  of  Attack/  Comparison  belZetn 

theCSachSBoxSMe?hoSalCUlati°nS  by  TW°  Forms  of 
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without  account ingCfJJ*thersingul«rXdownthehP5f3SUres  found 

!i^°r&  ex~ 

discrepancy  can  be  attributed  entir?iv  Jfach  boxes.  This 

has  to  satisfy  the  zero -pres sur£eo^d?ff  the  f*ct  thllt  on« 

(7)  ,  (8)  ,  rather  than  at  thf  j  tion  5*  P°ints  (sj  ,  ,6 , 

had  been  used  at  the  side  edges  with  cente^  h*lf,  M*ch  boxes 
if  boxes  were  half  boxes t  hhf  \LW17  .Centers  on  that  edge  fi  e 
proved.  However,  the  latter  a?^nefffn^C?uld  h*ve  beenta- 
from  the  practical  standpoint.  8em€nt  is  very  inconvenient 

calculated^^Procedur^fl^exhibit^6  ValiJ!:8  i2  T*ble  Iv-2 

vergence,  as'more  cho^dwise  boxts  lrt  faster  con- 

hy  the  characteristic -box  method.  **  taken'  th*n  those  found 

t?l„sBty„|1?Su0i?igeefo™US0|orqt31;e(5iil  *Vd  <42>-  «•  »•>- 

Pitchlng  moment  .bout  the  *  'mIs  n.2  ??lonle*!  Uft  "d 
half -span  (See  Fig.  IV. 1):  *  1  per  nit  sP*n  for  the  left 


n/  |  _  sfoai) 

*  WP-  - 


&"*»?/«*,/- S/k« 


=  - 1 


z. 

3jc 


+ 


fi 


“  S  £.£2?  correspondingSchordwise  tV*6 

catried  out  nun.rtc.llv 
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discrete  points.  In  general,  these  approximate  integrations 
if^nW06  *rrors  which  can  be  quite  significant,  especially 
f!W  P?ints  are  available  in  the  ranges  of  the  inte¬ 
grations.  For  instance,  at  Section  A-A  (Fig.  IV. 2)  the 
trapezoidal  rule  yields  />'??/>/) 


+  £f~' + *&+ 2 A*  +a9] + #4" 

Eqs.  (4 . I2a-b) 

At  Section  A-A  (Fig.  IV. 5)  the  rectangular  rule  yields  -/) 

^  -  j-  [-/ + h, 'i-  Az +fi,3  +h!-  +As ] 

^  %-A  =  &+{&¥)*&$]  Eqs.  (4.13a-b) 


- /  «necK  cneir  exact  counterparts.  Thei 

numerical  integrations  can  be  improved  appreciably  if  one 
accounts  for  the  slope  discontinuity  in  the  pressure  distribu- 

e£?sta£™S,-£he  }ln?  “  (Fi«-  IV  2)'  is  knSS  to 

the the  analytical  solution,  Eq.  (4.2).  The  effect  of 
the  infinite  slope  in  the  stream  direction  can  be  included 
as  shown  in  Appendix  C,  Ref.  31,  provided  the  fluctuations 
bhe  pressure  as  a  function  of  x  are  not  excessive.  This 
qualification  limits  the  refined  integration  technique  of  the 
above  reference  to  cases  where  the  leading  edge  is  supersonic, 
because  subsonic  leading  edges  produce  severe  fluctuations  in 
the  chordwise  pressure  distribution  as  estimated  by  any 
numerical  approach.  y  y 

For  the  particular  numerical  example  at  hand,  even  better 
numerical  integrations  can  be  devised.  For  example,  at  Sec¬ 
tion  A-A  (Fig.  IV. 2),  if  one  uses  the  trapezoidal  rule  for 
segment  ,  the  refinement  of  Ref.  31  for  the  sharp  drop 

in  pressure  in  segment  (3t)-(/z)  ,  Simpson's  rule  for  segment 

^ 9 )  bhe  trapezoidal  rule  for  segment  ( '9)  -s  one 

?ieldsS  satisfactory  accuracy.  This  combined  approximation 
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ft 

y  s  j  .  ^  J  .  -* 


£(b'  !£l) 

*"•9  8/ 


Tv.  1  8  E<ls'  (4.14a-b) 

preferabl^to^qs  fC*4  f°2) -(4°i3?qSa  definitely 

Jb^iled  h  «  «re  rational  evaluation  ^%Klnpfoved  E°™Us 
Cained  by  accounting  for  the  do™«sh  stelu&ar*"' 

In  Figs.  IV.6-IV  7  X' 

about  the  leading  *a*L  *  t"e  sP*nwise  running  lift 

especially  true  of  the  ehJSf  5*°?  considerably.  This 

IB^tS  satis^action  of  the  zern*’61'  Sti,C  ^rid  system,  which 
edges.  he  zero-pressure  condition  on  theside" 

o“J  explicitly,  ^ntt  i°™  that  the 

dovn™sh  and  for^'ch  o^rt.ricje?^ tic"?18'  v«1«i°"bin 
Jll%b  Slzes  *re  directly  relftpH1Jtifubox  s7stems.  where 
contrast,  if  one  employs  thf.  *  d  to  the  Mach  number  Rv 

U?Cdiff?ef?rehand^heyMach  number6  8ffd  ZyStem'  one  ^st* 

thl  is  encountered  to  "°  Par*icu- 

UritTfp'  ^ l«eJe^tsth?oei^UrthSu5e 

phragmnr^fr*  i^*^  °r”h^0*®^W*^°nErol^po?^  emanating  5  from' 
pnragm  region  in  an  irrponiiv  i*-ontroi  points)  cut  the  dia- 

planforms  with  subsonic  edjn.8rid  SysteiD  cannot  be  used  for 

steady  flow'at*  AiiCf “ngular  wing  at  an  angle  of  attack  ot  i 

,b“  (\*> . 13  1"fl'-«ced  Jy  bo?”  ("  “f-  *J-?-  ‘  dlsphraga, 
Y't-A  ls  influenced  by  (’faW/',  \*nd  (',«■).  In  tSrn 
gently  if  ehe  pressures  ere'to  be'^d.*;  <lnd  30  »”•  Cons™: 

U  b<Mte8  th*  *"•.  «"  infinite^nurtber  ?t fi^TSS* 
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Fig.  IV. 7  Spanvlse  Moment  Distribution  on  s  Rectangular  Wing  in  Steady 
Supersonic  Flight  at  Constant  Angle  of  Attack,  as  Calculated 
by  the  Methods  Indicated  on  the  Various  Curves 


CONFHKNTML 


Fig.  IV. 8  Square  Grid  System  for  a  Rectangular 
Wing  in  Steady  Motion  at  Constant 
Angle  of  Attack.  (  A J  =  /.  z  ) 


In  d;*,  d  do^n«a8he8  on  these  boxes  determined, 

in  practice,  the  number  of  boxes  must  be  limited,  so  that 

;ii^wibey°nd  8T®  b°undary  line  off  the  wing  must  be 
arbitrarily  assumed  to  have  zero  downwash.  As  a  start,  let 

the  following  boxes  and  those  which  are  further  out  have 
zero  downwash 


The  nonzero  downwashes  etc.  are  then  determined 

by  the  conditions  of  zero  pressure  at  the  centers  of  boxes 
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’  •  ,Here  *  sequential  solution  is  not  pos- 
ble,  since  laterally  adjacent  boxes  affect  each  other.  The 

fhaJ  yfel?  th®  doynw**he*  for  the  diaphragm  boxes 

J!  Jked  •im4ltaneou«ly-  The  resultant  pressures  at  three 
points  on  the  planform  are  computed  to  be 

~°- 573  H  '  £>'2=  ~a-s 73  %,  4'  =  -  56  , 

as  compared  with  the  exact  values 

/>,  -  ~fl.  605,92,  -^.605  92,  =  -a.3  n,GS  . 

IJurin8  this  computation  the  condition  of  zero  pressure  on  boxes 
gives  *  (Z'c)  ’  eCc*»  violated.  In  fact,  a  simple  computation 

ti*  =  C- ¥'9' ¥5'  k,c  =  *  oc 

Next,  one  additional  box  is  taken  in  each  row,  i.e.. 


4C 


o  ' 


j  ^  o  . 


The  corresponding  pressures  at  two  points  on  the  wing  turn  out 
to  be 


h,  /  m  -  *°7,  5b,  -  -  0.V07,S& 

If  the  foregoing  process  is  continued,  taking  one  additional 
box  at  a  time  in  each  row,  the  results  can  be  shown  to  deviate 

nK>5VnVBore  fro“  the  e**ct  v*lue«-  Therefore,  the  first 

”iere  only  diaphragm  boxes  with  centers  inside  the  Mach 
line  ab  are  considered,  appears  to  give  the  best  possible  ac¬ 
curacy,  even  though  the  pressure  conditions  at  (/S)  ,  />  c) 
etc.,  are  violated.  7  ’ (  '  y  ’ 


„  .  A ,81jP>r  difficulty  occurs  for  planforms  with  swept  super- 

^3eadln8  6<lfe8  ^c,£  »  P°int*  *■  «nd  S-  ,  Fig.  III. 3  where 
according  to  the  box  scheme).  While  this  violation 
or  the  zero-pressure  condition  is  not  serious  for  supersonic 
leading  edges,  as  indicated  by  the  results  of  Table  I II. 3  it 
leads  to  much  larger  errors  when  there  are  subsonic  edges.  As 
cen  be  seen  from  Table  IV. 3,  the  pressures  for  the  first  case 
treated  above  diverge  more  from  the  exact  for  boxes  farther 
downstream  and  fluctuate  violently. 
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Procedure  (1) 
Using  Square 
Boxes 


pi,l  -0.573,96  -0.605,94 

P2, 1  -0.262,56  -0.311,66 

p2, 2  -0.573,96  -0.605,94 

P3, 1  -0.310,63  -0.237,34 

p3, 2  -0.356,47  -0.288,28 

p3,3  -0.644,50  -0.605,94 

p3,4  -0.808,53  -O.827.89 

p4 , 1  -2.410,91  -0.199,21 

p5,2  -2.652,43  -0.311,66 

_ p6,3 _  "0-098,30  |  -0.370,07 

Table  IV. 3  Pressures  on  the  Mixed  Region  of  a  Rectangular 

At?Lw?teady  ?°ti0n  at  a  Constant  Angle8offlr 

Attack  (Square  Boxes,  /*)=/.&  ) 
fhe  Side -Edge  Problem  in  Unsteady  Flow 


Exact  [Eq .  (4.1)] 


p3, 1 


reUttvely£“^le^robl“f  h?Ye  with  th< 

the  general  P^cedGres  jnsfoutiined  «e '  y  T0*1? 

dition  to  regular  ATP  p.kuo  u  are  still  valid.  In  ad- 

/  p\  of  Se:fSJA?5  VIC  ta°les’  however,  the  use  of  Procedure 

wash^tstrfbutior^lrjL'nLvlii1*" 

for  the  pressure  at  the  point  (  r Vo  J  i.o’VinV'L0^8 
wa3h  functions  given  by  Eqs.  (4.7a-b)  UC  t0  the  down" 


P3,4 
p4 , 1 


P6,3 


-0.310,63 

-0.356,47 

-0.644,50 

-0.808,53 

-2.410,91 

-2.652,43 

+-0.098,30 
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(origin  on 
Eqs . 
(origin  on 


diaphragm) 

(4.15a-b) 

planform) 


Here  is  the  amplitude  of  the  singular  downwash  strength, 

now  assumed  to  have  simple  harmonic  time  dependence,  co  is  the 
circular  frequency  of  oscillation,  and 


r 


J  v^/s-sj  yps 
0  L  / 


(*-V)  ■ 


Eqs.  (4.l6a-b) 


The  notation  is  that  of  Fig.  IV. 4. 

By  suitable  changes  of  variables,  the  double  integrals  de¬ 
fining  jr  and  J*  can  be  reduced  to  the  single  integrated  series 


Eq.  (4.17a) 
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■Pn  =PJ*i<rJM)=e 


cea 


in 
£  £ 


vT 


{iKrf-Z'fr-'r  )I[e  [**--1 fj] 


where 


,  r,  a>M 

/  up 

-=f 


Eq.  (4.17b) 


SinceT  iS  ^  BeSSel  function  °£  the  first  kind  end  order  /> 

Jo(o)  -/  J  (o)  —  <9 

Eqs.  (4  I7a-b)  in  conjunction  with  Eqs.  (4.l5a-b)  confirm  the 
expressions  previously  given  for  steady  flow  [Eqs.  (4.8a-b)]. 
The  infinite  series  of  products  of  Bessel  functions  in  the 
integrands  are  rapidly  convergent  for  the  frequency  and  Mach 
number  ranges  of  interest,  so  that  only  a  few  terms  need  be 
retained.  The  integrations  indicated  by  Eqs.  (4.l7a-b)  must 
be  carried  out  numerically.  In  this  connection,  it  is  worth 
noting  that  the  integrands  are  functions  of  z*  only.  More- 
over,  for  the  Mach  and  characteristic  grid  systems  (the  only 
ones  which  permit  introducing  the  singularity  effect),  rr  is 
an  integer  or  a  rational  fraction,  and  5=  is  a  rational  frac¬ 
tion  less  than  unity.  A  suggested  method  for  the  evaluation 
of  these  integrals  is  given  in  Appendix  A. 


.  .  an  exe?nPle  to  demonstrate  that  the  inclusion  of  the 
side-edge  downwash  singularity  is  feasible  and  results  in 
fairly  accurate  airload  distributions,  calculations  have  been 
carried  out  for  a  rectangular  wing  in  rigid  body  translational 
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oscillation.  Since  the  downwash  is  constant  everywhere  rh*. 

orjM^"1!tlVrt'm  C*n  b*  “««•  without  «count^rihe 
problems  caused  by  -^“-variation  over  a  box  This  svst-em 

cSndi?lnlnCf  U  peraits  the  satisfaction  of ' the^en/pressilre 
condition  at  points  along  the  side  edge.  In  Figs  IV  9-IV  lo 
the  aerodynamic  derivatives  (Ref.  13)  8  9 IVl°- 

.  S/>a» 

L  +  c  L,-  - - - - - - 

M  V  L  M*  =  about  eJje sikt/) 

arassg. 

cop1?  S". •  V V^°{s  •>«,"b1^1??e"^asya?|?^e«- lth 

the  amplitude  ^  of  the  oscil lifinn  tk^  . 

amplitude  is  cne^ oscillation.  The  constant  downwash 

it?  —  cu>$B 

13)  thShi<>h'ne*T  £”«b"s*of  Nelsona°Raineynand^WatkinsS(  Kef’’"1*'1 

■SSSrSSS&SSsaescS- 

Z-  -*4*1.  1.913, 

R?fiSnU°¥Sl|(Under  bh!  fropuoncy- 
Ref  «  -b°d  8ives°values  %  ^ 

two-dimensional  region  which  are  a  few  per  cent  off  fr™  eu. 
(RlfeS36)'dinfhe?tforealUeSnCCOrdin8  to  Garrick  a"d  Rubinow 

h£~*  Stt,”  »™"ssk  Exrx  sag. 

;s.  •?«:  ss.*\ri:  sya  ras-s&s,. 
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TN  3076  SERIES  EXPANSION 


- EXACT  —  TWO  DIMENSIONAL 

O  4  CHARACTERISTIC  BOXES-WITH  SINGULARITY 
•  4  CHARACTERISTIC  BOXES  -  NO  SINGULARITY 


Fig.  IV. 10 


WADC  TR  56-97,  Part  1 


?i!TS?ior‘}e8£  SPa?wi?e  Moment  Distribution  on  a  Rectangu 
);?r  in  translational  Oscillation,  as  Calculated  by 

the  Methods  Indicated  on  the  Curves.  (  &  *  o.  3  ^  =J.Z  )  Y 
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11D  .  In  5he  calcul-ation  above,  only  four  chordwise  boxes  were 
used,  yet  accuracy  sufficient  for  most  purposes  was  obtained 
By  contrast,  when  the  same  number  of  boxes  were  used  and  the* 

p?”lslon3reluuIdty  h‘d  bee"  *CCOUnte<1  ™<=»>  Poorer 

the  pressure^wes^checke^only  ll™ 

one  point.  For  point  14,  Fig.  IV. 2,  °  Ay  aC 


4L 

p 


=  —  O.  307, &0  -t-o./lb,  Wi 
^  -  0 .  312, 60  i-0. 1/5J  /&  c 


( faf.  13) 

^  (TAarac  7tr/s//c.  ^oxes 

u,,#,  ,«&»»****). 


the  technique  °!  including  the  downwash  singular- 
ity  was  shown  to  be  feasible  for  unsteady  flow,  some  question 
arises  as  to  its  practical  value.  In  steady  flow  theAIC's  as¬ 
sociated  with  the  singular  distributions  have  such  simple  mathe- 

S^reC3}ff?P.r?^SihnS  lEqS‘  -8«-b)  ]  that  Procedure  (2)  is  no 
more  difficult  than  Procedure  (1).  When  the  motion  is  unsteady 
however,  one  must  first  tabulate  the  functions  ln  and  ^  [Eqs/’ 

(4 . I7a-b) ] ,  each  of  which  depends  on  three  variables.  The 
first  g lane e*1  for * t he  ^  f o 1  lowing^eas  ons  T***"  “  “  —  at 

(1)  Tables*  need  be  prepared  only  for  discrete  values  of 
<r  and  cr  .  When  a  large  number  of  chordwise  boxes 
s  taken,  the  ranges  of  cr  and  cr  may  become  exces¬ 
sive.  However,  it  will  be  shown  below  that  the  in¬ 
clusion  Oi.  the  downwash  singularity  is  no  longer 
necessary  when  a  great  many  boxes  are  used. 


(  2 )  1n  and 


A 


are  both  smoothly  varying  functions  of 

the  frequency  parameter  €  .  Therefore,  it  is  not 

necessary  to  take  a  large  number  of  ^  -entries  be- 

intervalserP°^ati0n  C<m  ^  done  accurately  over  large 

Since  tables  are  not  as  yet  in  existence,  it  is  desirable 
to  establish  a  lower  limit  on  the  number  of  chordwise  boxes  re- 

?S15£hnS^Ih?  Slde  edgeKif  Pr°cedure  C1}  i*  to  be  employed 
If  Mach  boxes  (or  square  boxes  for  are  used,  it  can  be 

^own.Lc.f.^Figs.  IV . 6 - IV . 7 )  that  the  significant  ina^c^acies 

•These  tables  are  three-dimensional,  wherea*  the  tabulations 
for  the  usual  AIC  are  four-dimensional. 
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in  the  lift  and  moment  distributions  are  confined  to  the  first 
two  or  three  stations  near  the  side  edge.  The  more  boxes  that 
are  taken,  the  more  these  critical  stations  will  be  confined 
to  a  small  portion  of  the  total  mixed  region,  and  the  lifts  and 
moments  will  then  be  closer  to  their  exact  values.  This  is 
demonstrated  by  the  points  calculated  by  the  ten-box  analysis 
in  Figs.  IV.6-IV.7-  There  the  total  lifts  and  moments  for  the 
mixed  region  are  expected  to  be  sufficiently  accurate.  On  the 
basis  of  this  steady-state  example,  a  minimum  of  ten  chordwise 
boxes  is  believed  to  be  required  when  using  Procedure  (1). 

To  substantiate  this  conclusion  for  unsteady  motion, 
another  numerical  example  has  been  carried  out.  The  lift  and 
moment  distributions  calculated  with  six  and  ten  square  boxes 
along  the  chord  for  a  rectangular  wing  in  rigid-body  pitching 
oscillation  about  its  mid-chord  ,  M=/.5  )  are  presented 

in  Figs.  IV.11-IV.12  and  compared  with  corresponding  three- 
dimensional  values  from  Ref.  13  and  two-dimensional  values 
from  Ref.  36.  For  this  combination  of  Mach  number  and  reduced 
frequency,  the  former  can  be  regarded  as  a  standard  of  refer¬ 
ence. 


It  is  seen  that  with  ten  square  boxes,  the  AIC  results 
are  much  improved  over  those  obtained  with  six  boxes.  The 
total  lifts  and  moments  over  tbe  mixed  region  with  ten  boxes 
are  acceptable,  since  the  major  discrepancies  are  confined  to 
a  limited  area  near  the  tip.  If  the  contributions  of  the 
mixed  regions  to  the  airloads  are  small  compared  with  those 
of  the  total  wing,  one  might  tolerate  somewhat  poorer  accuracy 
in  these  regions  and  accept  fewer  than  ten  chordwise  boxes. 
However,  this  can  be  done  only  for  certain  special  planforms 
and  motions.  If  the  motion  involves  spanwise  elastic  deforma¬ 
tions  such  that  the  deflections  near  the  tip  are  large,  which 
is  often  the  case  in  flutter,  the  aerodynamic  loads  will  be 
concentrated  on  the  outboard  sections.  Maximum  accuracy  is 
then  demanded  in  the  mixed  region,  regardless  of  the  aspect 
ratio.  Therefore,  the  conservative  and  generally  acceptable 
rule  seems  to  be  either  to  adopt  at  least  ten  boxes  along  the 
wing- tip  chord  or  make  special  provisions  to  include  the 
proper  downwash  singularity  off  this  tip. 


IV . 3  The  Subsonic  Leading  and  Trailing  Edge  Problems 

In  the  previous  section  a  method  was  devised  which  in¬ 
cluded  the  effect  of  the  downwash  singularity  at  a  side  edge  on 
the  pressure  distribution  over  a  finite  wing  oscillating  in 
supersonic  flow.  When  a  similar  procedure  is  attempted  for  a 
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- EXACT  TWO-DIMENSIONAL  TR  846 

- —  SERIES  EXPANSION  TN  3076 

X  SIX  SQUARE  BOXES  (PROCEDURE  I) 
O  TEN  SQUARE  BOXES  (  PROCEDURE  l) 


Fig.  IV. 12  Dimensionless  Spanwise  Distribution  of  Pitching  Moment 
About  the  Mid-Chord  of  a  Rectangular  Wing  in  Pitching 
Oscillation  About  Its  Mid-Chord  Axis,  as  Calculated  by 
the  Indicated  Methods.  m  =  /.5  ) 
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.■Scb0cnoSpu!Jt^„1.Chtb1.%n?ht..P*:'»Jel  iV*  fl0“'  “  le.d.  to 
if  not  impossible.  *  PProach  i«  regarded  as  impractical 

edge  ~.”%P?e1^ee,deStetd1bnvS.eCbtr1°;.IIf(1  *?“  *  —P*  leading 
u la ting  the  AIC^s  for  partial  boxea  1if*th?.°rd*r  *?  *vold  tab' 
accepted,  the  question  arises  *!  thi8  aPPr°*i®ation  is 

lar  downwash  distribution*  uh  If  one  can  represent  singu- 

the  axis  of  the  singularity  parallels ^hf*.  °?  *5®  actuaI  win8 

form  with  the  jagged  edges^!!  ko  th  ept  edge-  To  con- 
adjusted  .o  thit*ft  foil".  ?he  rhi“  *xl<  »«  be 

line.  Practical  difficulties  are  segments  of  the  broken 

assurance  that  the  distortion  of®i-knf:oynteTed'  and  there  is  no 
improve  the  results  over  £hat  Soul^h! ijguiar  tiistribution  will 
singularity  altogether  "  *  ld  be  found  br  neglecting  the 

scheme  for88ubsonic8edgeshotherUthabilJHy  th®  singularity 
vestigated  another  alternati™  2l? u®ide  edges,  the  authors  in- 
aingularity  by  a^llticll  Conaiata  of  removing  the 

tion  of  the  solution  by  the  box  method6**???*^6  re®ainin*  Por- 

Tbe  downwash  singularity  at  a  noJnrr°V®M*®qUaily  im* 
leading  edge  of  a  finite  wing  is  deter^iiLd  the  subsonic 

portion  of  the  wing  inside  the  forwiJd^taS!  th®  Motion  of  the 
this  point.  The  influencing  Irea^Jv^!?? h  lines  emanating  from 

ing  edge,  if  the  latter  is  also  auSLJj^1U2®  *  P*rt  of  the  trail - 
method  for  calculating  this  Bineulllitt  H*nce>  an*  analytical 
just  as  well  to  find  the entirl^loId^Li n  advance  could  be  used 
aede  the  AIC  approach  entirely.  *  distrlbution  and  would  super- 

the  conclusion  tbst,  if 

subsection.  YP  represented  b7  Procedure  (2)  of  the  previous 

pl  informs  that  progressively^ncreasing^he156  Ph°V®dffvr  other 
the  chord  causes  the  result?  for  the^fne^lized6^  b°?ea  along 
verge  toward  the  exact  values  To  I.SIkwS.1**?  forces  to  con- 
number  of  boxes  required  to  ichilSe  J^tiiff^®  l0W®r  J1?1*  on  the 
treatment  of  a  subsonic  leadinvadL aatiafactory  precision  in  the 
merical  example  is  first  oregpJfprf8®  problem,  fbe  following  nu- 
with  sweep  angle  A  y.P/  ^  triangular  or  delta  wing 

sweep  angle  A-  cd  '**<  is  considered.  The  flow  is  at 

of  et  t.(£  j:  ( TcthlAinhi‘  '^:dLT!tlon  «  «—t-t  angle 
retical  airloads  can  be  derived  *_For„5bi8  caae»  tbe  exact  theo- 
vti«ic. ........ i  glv,„ ..  th* 
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cotr  -A-  j 


L  ft /uru7  sfean  _  Cot  A-  /  /  _/tt  f 

(zsu  °c)2^  _  ~ 


Moment /u„,t  s/xw 

7*fv5j&ir~= 


cot  A 


Eqs.  (4.l8a-c) 

?ure  ^  if  the  comPlete  elliptic  integral  of  the  second  kind, 

i?J,.£?£dln*teJ  4K  *Ze  thoSe  of  FiS*  Iv*  13»  and  the 

pitching  moment  is  about  the  y  -axis:  cotA/e=  =  W  9^ 

By  spanwise  integrations  of  Eqs .  (4.l6a-c),  one" obtains  for  the 
dimensionless  total  lift  and  total  moment 


/_  _ 


■JUCot^L 


-0.987,09 


Moment  JCcot  z/i 

(2fl/tt)[ZS3~)  3  ~E  = 


/•  3/4 II 


Eqs.  (4 . 19a-b) 


lue  filure  illustrates  the  pattern  of  Mach  boxes  adopted  for 
the  AIC  calculation.  At  this  particular  Mach  number,  inciden¬ 
tally,  Mach  and  square  boxes  are  identical. 

In  each  box  of  Fig.  IV. 13  two  values  for  the  pressure  are 
given,  the  upper  being  that  calculated  by  the  box  scheme  and 

1^er,f5?n’^he,exact  forinul*  £E<1-  (4.18*)] .  It  is  seen 
that  the  individual  pressures  given  by  the  numerical  method 
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Fig.  IV. 13  Dimensionless  Pressure  Distribution  on  a  Tri¬ 
angular  Wing  at  Constant  Angle  of  Attack. 

\A  =  0/  (Lot  A.  ~0.  6^  s/e) 
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£hU«tU*te  conaider«bly  in  the  stream  direction.  A  similar 
phenomenon  was  observed  in  connection  with  a  sweet  suoerlLf^ 
edge  (see  Section  III.l),  although  there  it  w«  S£ch  le« 
severe  Despite  these  fluctuations,  the  spanwisriife  !£d 
rSSS*  ^*«^«lon.  follov  th.  treid. 

?ho:.1Xf„rTS«r^;.for  th'  Uft  *PP'*r  •— it  better  than 
seem  SSl  UfC  “J  Cot<1  "»«« 

^  T!?e>,-0nVer8?nceui8  not  unifor»  but  proceeds  in  cycles 
It  happens  in  this  example  that  the  error  increases 

*nd  then  ,uddenly  jumps  to  a  smaller  value 
h™.PI,en?“en0n  r*P®atinS  each  time  the  number  of  ' 

Sf  Hri  increased  by  three.  The  only  rational  way 

!  conver8ence  is  to  examine  the  maxing 
deviations  of  successive  cycles. 

rk*  convergence  based  on  successive  cycles  is  slow 
Iierf£ore:  ®*ny  ®ore  than  twelve  boxes  would  have  to 
S-  t0  achiev®  a  significant  reduction  of  the 
of  Table"  IV™0**  characteri8tic  of  ^e  lower  portion 

S£?StiCf1  Con!iderations  of  manageable  computations 
thi  i^kC  any  increa?e  far  beyond  twelve  boxes,  since 
%h?«°r  8^OW®  *3  1?88t  in  ProP°rtion  to  the  square 
of  this  number.  It  is  therefore  tentatively  con- 

£orit  .I**  tWalu f  boxes  alon*  the  “idspan  chord  af- 
acceptable  compromise  between  accuracy  and 
computational  labor.  Two  additional  examples  of 
triangular  wings  with  subsonic  leading  edges  in  un¬ 
steady  motion  are  presented  in  SectioS  V.8  It  is  seen 


*In  all  these  cases,  the  box  sizes  are  adiusted  so  th*i-  i-ho 
^^n,diite*a*^**e^a^h*oC^b®1*^®*b°fe*h®^wingr€PSincetthe^>^an” 

..:s  te'.K'ssui'Ss:;;,!!:."—1"  “• 
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Number  of 

Chordwise 

Boxes  at 
the  Midspan 

Station 

L 

'  1 

-1.28 

-1.536 

2 

-1.28 

-1.792 

3 

-1.028,76 

-1.193,55 

4 

-1.075,26 

-1.378,69 

5 

-1.130,02 

-1.540,59 

6 

-1.043,96 

-1.345,13 

7 

-1.061,81 

-1.399,60 

8 

-1.093,96 

-1.483,04 

9 

-1.040,12 

-1.365,48 

10 

-1.048,12 

-1.389,83 

11 

-1.069,35 

-1.443,09 

12 

-1.035,23 

-1.369,79 

Exact 

[Eqs.  ( 4 . 19a-b) ] 

-0.987,09 

-1.316,11 

i.abie  IV . 4 .  Total  Lifts  and  Moments  on  a  Triangular  Wing 
Constant  Angle  of  Attack  in  Steady  Flow,  Cal¬ 
culated  by  the  AIC  Method  with  an  Increasing 
Number  of  Chordwise  Boxes.  M=Vsl) 
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there  that  six  or  eight  boxes  along  the  roidspan  chord, 
which  proved  sufficient  for  purely  supersonic  plan- 
forms,  are  inadequate  when  subsonic  edges  are  present. 

If  large-scale  digital  facilities  are  available  and 
costs  do  not  have  to  be  held  down,  the  requirement  of 
higher  accuracy  in  particular  cases  may  increase  the 
tentative  minimum  by  a  substantial  amount.  The  ex¬ 
perience  of  the  authors  to  date,  however,  does  not 
permit  any  quantitative  estimation  of  the  probable 
error  associated  with  any  given  number  of  boxes. 

Originally  it  was  believed  that  one  could  alleviate  the 
errors*  in  calculations  like  the  foregoing  (wings  with  super¬ 
sonic  trailing  edges  and  subsonic  leading  edges)  by  making 
use  of  the  reverse-flow  theorem.  After  reversing  the  flow, 
it  is  a  case  of  dealing  with  a  planform  with  subsonic  trail¬ 
ing  edges  and  supersonic  leading  edges,  and  thus  eliminating 
regions  of  pressure  singularity.  It  will  be  demonstrated  that 
this  artifice  is  unsuccessful.  The  identities  implied  by  the 
reverse-flow  theorem  (Ref.  17)  for  the  total  lift  and  total 
moment  on  a  wing  at  constant  angle  of  attack  in  steady  flow 
can  be  written 

L=JJb-' 

M=  2f  JJ/> =  ^ JJh  *  EqS  '  •  20a’b) 

4-  £ 

Here  and  $>_  are  the  regions  of  integration  over  the  plan- 

form  in  forward  and  reversed  flows,  respectively;  is  the 

pressure  due  to  a  downwash  L**.  in  forward  flow,  is  due  to 

a  downwash  bU  and  is  due  to  a  downwash  (*-/2.6-)~  1/  in  reverse 

flow.  If  the  assumption  of  constant  downwash  over  each  dia¬ 
phragm  box  near  a  subsonic  trailing  edge  were  less  critical 
than  the  same  assumption  near  a  subsonic  leading  edge,  then  the 
reverse -flow  approach  would  improve  the  accuracies  of  the  total 
generalized  forces.  However,  if  one  considers  the  triangular 
wing  and  solves  the  same  problem  in  forward  and  reversed  flows, 
one_obtains  identical  results.*  This  indicates  that,  as  far  as 

*It  is  not  obvious  that  the  two  configurations  should  yield  the 
same  results,  even  though  Eqs .  (4.20a-b)  are  exact  expressions 
from  linearized  theory,  since  solution*  by  the  box  scheme  for  both 
flow  directions  are  approximations. 
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f°rfe  estln,ations  by  the  AIC  method  are  con- 
sub^oMr  ^fl^eren5  err°rS  fntroduc®d  in  the  treatment  of 
sonlc  UadJng  aSles  885  JUSt  aS  large  as  those  from  sub' 

_  .  i£  one  admits  that  a  minimum  of  eight  boxes  along  the 

is-re<luired  for  purely  supersonic  planforms  and 
twelve  boxes  for  planforms  with  subsonic  leading  edges,  it 
JSOnabie  to  extrapolate  and  set  a  minimum  of  sixteen 
in£hJiS£r  fwr  with  both  subsonic  leading  and  trail- 

bee  dgj  No  calculfltions  on  examples  of  this  sort  have  yet 
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SECTION  V 

SUPPLEMENTARY  NUMERICAL  EXAMPLES 

Certain  additional  examples  of  the  application  of  the  AIC 
method  to  oscillating  delta  wings  are  included  in  this  section 
Three  of  the  cases  involve  the  use  of  square  boxes  at  f>5  ; 
the  fourth  illustrates  the  Mach-box  scheme  at  /.Z  ,  which 
is  regarded  as  just  about  the  lower  practical  limit  of  linear¬ 
ized  supersonic  theory.  These  examples  were  worked  out  in  the 
course  of  investigating  various  difficulties  associated  with 
numerical  methods,  and  for  the  sake  of  clarity  were  left  out 
from  previous  Sections.  The  reader  will  recognize  from  ac¬ 
companying  discussion,  however,  that  each  result  tends  to  sub¬ 
stantiate  one  or  more  of  the  conclusions  and  rules  stated  in 
foregoing  sections. 

V. 1  Purely  Supersonic  Planforms 

Case  1:  The  spanwise  lift  and  pitching  moment  distribu¬ 
tions  are  given  in  Figs.  V.1-V.2  for  a  triangular  wing  with 
leading  edge  swept  back  forty-five  degrees  flying  steadily  at 
and  constant  angle  of  attack.  Three,  six  and  eight 
square  boxes  were  employed  in  the  computations.  All  these 
airload  distributions  are  remarkably  close  to  the  exact  ones, 
which  also  appear  on  the  figures.  One  obtains  even  more  ac¬ 
curate  estimates  for  the  total  lift  and  total  pitching  moment 
because  of  the  averaging  effect  of  the  spanwise  integrations. 
Although  in  this  particular  case  few  boxes  yield  excellent 
results,  the  same  cannot  be  said  for  other  purely  supersonic 
p/.anforms  undergoing  elastic  deformations  at  high  frequencies, 
as  shown  in  Section  III. 

Case  2:  The  same  wing  as  Case  (1)  is  considered  to  per¬ 
form  parabolic  bending _oscillations  at  ^  =  o.Sl,  this  motion 
being  described  by  £■=  (jf/ztfe. tco  ,  where  is  the  tip 

amplitude.  With  eight  chordwise  boxes  at  the  root,  one  obtain 
excellent  estimates  for  the  total  lift  and  pitching  moment 
about  the  apex,  as  shown  by  the  following  comparison  with  the 
corresponding  quantities  from  the  exact  theory  of  Ref.  (17): 
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Attack  in  Steady  Flow.  (A  =  n-s 


total  moment 
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comum 


Total  Lift: 


Total  Momenr 


=  -0.023,09  -  0.158,831 
-  -0.023,21  -  0.158,911 

*  -0.038,06  -  0.252,951 

*  -0.038,17  -  0.251,201 


Exact 


Box  Method 


Exact 


Box  Method 


substentlatlo^h|0£thc8quallt«lS  "ay1be  “Id  to  add  further 
be  taken  along  the  «idSp.“£^%Fj%fJ«  *1*1*  hoxe.  should 
V . 2  ,u-a— .  . .  ....  .  S  ^  suPer9°nlc  edges 


Planforms  with  Subsonic  m,.-. 


are  given  In  FlgsnWVS|-v16tfo?d.PJtfhlng,moment  distributions 

adge  swept  back8slxty  degrees  perfo™?^Uf  "ln«  wlth  la«d?n| 
tion  anti  pitching  oscillationsPabmi^h8  rigid_body  transla-8 
*t  &-O.Z,  and  M  =  /.s  .  giebf-  ;.Ut  root  mid-chord  axis 
along  the  midspan  chordiine  8  it  9“are  boxes  were  employed 
discrepancia3  exist  betieen' 'the  ai-lSSJT  thaf  aP^clable 
and  the  more  exact  results  of  R.f  !^d  froin  the  AIC  scheme 
moment  distributions.  Although  5k  20 '  esPec*ally  for  the 

ofrtheUt  t0  Somewbat  better  becluse^f  thftS  *nd  ®on,ent:s 

SSl^SoSf  if SteT^'t^^g 

sufficient  increase  in  precision1?118  the  midsPan  chord  yields  8 
putational  labor.  ThLS  e  the  additional  corn- 

even  more  significant  for olimSSS  J  be  ?xPected  to  become 
tions.  Further  enlarging  the  nSlr  Je!8°ing  elastic  defonJsa- 
oes  not  seem  to  be  justifiable,  hov?everb°XeS  beyond  twelve 
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Case  4:  The  spanwise  lift  and  moment  distributions  are 
given  in  Figs.  V.7-V.8  for  a  triangular  wing  with  leading  edge 
swept  back  forty-five  degrees,  at  ■6=0.1,  and  /Y=/*£,  executing 
a  spanwise  parabolic  bending  oscillation.  Six  chordwise  Mach 
boxes  were  employed  in  the  calculation.  The  accuracies  ob¬ 
tained  for  these  airload  distributions  are  seen  to  be  rather 
poor.  As  in  the  previous  case,  one  must  apparently  resort  to 
the  order  of  twelve  chordwise  boxes  at  the  root  if  one  is  to 
get  reasonable  approximations  for  the  airloads. 
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SECTION  VI 

ALTERNATIVE  NUMERICAL  APPROACHES 

The  numerical  method  which  forms  the  principal  topic  of 
the  present  report  is  based  on  the  concept  of  the  direct  aero- 
dynamic  pressure  influence  coefficient,  which  gives  pressure 
at  a  point  due  to  an  influencing  elementary  area  with  a  pre¬ 
scribed  downwash.  Two  other  types  of  influence  coefficients 
have  been  proposed  and  were  studied  to  some  extent  by  M.I.T.: 

(1)  The  inverse  pressure  influence  coefficient,  which 
gives  the  downwash  at  «  point  due  to  a  prescribed 
pressure  distribution  over  an  influencing  elementary 
area  This  has  the  very  desirable  feature  that  it 
avoids  consideration  of  the  diaphragm  region,  where 
the  pressure  is  everywhere  zero,  in  cases  of  wings 
with  subsonic  edges.  It  forms  the  basis  of  numeri¬ 
cal  methods  in  unsteady,  three-dimensional  subsonic 
flow  (Refs.  28,  41),  where  the  direct  AIC  cannot 
readily  be  calculated. 

(2)  The  lift  influence  coefficient,  which  represents 
the  total  lift  acting  on  an  elementary  area  due  to 
a  prescribed  downwash  over  another  elementary  area. 
This  coefficient  has  been  found  to  lead  to  complica¬ 
tions  ( e .  g . ,  two  additional  integrations  must  be 
performed)  and  affords  no  special  advantage  over 
the  originally  proposed  AIC.  Therefore,  for  the 
sake  of  brevity,  numerical  investigations  based  on 
this  type  of  influence  coefficient  are  not  reported 
here . 

It  will  now  be  demonstrated  that  the  inverse  pressure 
coefficient,  which  seems  quite  promising  at  first,  can  be  used 
with  no  difficulty  for  steady-state  analyses.  For  unsteady 
motion,  however,  the  required  tabulations  become  so  complex 
as  to  render  this  method  quite  impractical. 

VI-!  The  Inverse  Pressure  Influence  Coefficient 

.  x2%*  recent  publication  Voss  and  the  present  authors 

(Ref.  28)  outline  a  procedure  for  determining  the  airloads  on 
an  oscillating  elastic  wing  in  subsonic  flow  at  low  values  of 
the  frequency  by  the  method  of  the  inverse  pressure  influence 
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sures  at  the  same  set  of  nn?n?.  tne  wing  m  terms  of  the  pres- 
usually  too™  in  ^advene  e  ^t  he d°esu  1  ting  EJsrf 

simultaneous  equations  Mat  ba  solved  for  the  S0«  pref- 
s«J  becw.i  £o?  ^*"8f  eht>  lndlr««  process  lsMces- 

SSiira?  boo£w,tie%1”tthef:,lu.‘r?±n  svr  &a>“  *> 

tion  gives  the  pres^redljectlvjj  1*  ln“*r«}  equa- 

of  the  weighted  downwash  [as  inEq  (!”))*  surf,ce  integral 


£*•  (6.1) 


verted ^o'yield  E’u*tton  <6-l>  «n  be  in¬ 
terns  of  a  surface  SEUJiVrt*  rfp"*'nt8  the  downwash  in 
tion,  C*  integral  of  the  weighted  pressure  distribu- 


4T 


Eq.  (6.2) 


where  /f  is  the  so-called  "inverse''  kernel. 

representatlor^oi^isturbed^reas^nenr'sub"*1!  t'*eH<'la*>'>r*8n' 

be  considered  ainei  fn  ?  near  subsonic  edges  need  not 

tally  zero  Mor««r  tS™  ff8fons  *5e  Pressures  are  identi- 
downwashes  at  subsonic  edges.  need  to  introduce  singular 


*(*-*'*-7)- -jj-Q  U(*-t-/>ini)^jqFftf  Eq 


(6.3) 


where  U  is  the  step  function 
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U(x)  .=  / 


-0  X~° 


Eq.  (6.4) 


If  the  inverse  pressure  influence  coefficient  is  defined 
as  the  dovnvash  at  a  point  due  to  a  unit  constant  pressure  over 
an  influencing  box,  then  in  the  Mach  grid  system  (Fig.  II. 3)  the 
relation  corresponding  to  Eq.  (2.14)  for  the  "direct  case"  is 


Fq.  (6.5) 


Using  Eqs.  (6. 2) -(6. 3),  one  can  derive  the  following  expres- 
sions  for  these  coefficients: 


P  =  -1  =  o 

0,0  ' 


-1/  > /t  -2*0 


Eqs.  (6.6a-c) 
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A  partial  list  of  numerical  values  of  is  given  in  Table 

J?*1*  is^ significant  that  this  table-covers  all  values  of 
inEqs  sinC*  this  Pararoeter  does  not  appear  directly 


0 

1 

2 

3 

4 

5 

0 

-1 

1 

-1.016,98 

0.508,49 

2 

-1.230,00 

0.485,75 

0.129,25 

3 

-1.254,91 

0.444,66 

0.117,67 

0.065,13 

4 

-1.263,09 

0.435,10 

0.097,68 

0.057,88 

0.040,89 

5 

-1.266,79 

0.431,07 

0.092,27 

0.045,53 

0.035,82 

0.028,70 

Taole  VI. 1  A  Short  Table  for  Inverse  Pressure  Influence 
Coefficients  in  the  Mach  Grid  System. 

Steady  State  (^=<5  )  7 


tnno.  R?§ardl«ss  of  which  box  scheme  is  to  be  followed,  it  is 

Thli  !^Ci?n?.t£>  !8SUxn?  f  constant  pressure  over  each  box. 
This  step  limits  the  tabulations  to  a  single  tvoe  of  inverse 

?hree8:SLeui^1T^e,COeffiCie?t*  0ne  ■ttst  carefully  examine 
y  of.  this  fssunption,  especially  near  subsonic 

iS?-ln8  edges  wh«r«  the  pressures  have  singular  behavior.  For 
coineS?^8*^  nun^rical  examples,  which  were  considered  in 

SSSTbSK  SbilSS^Stet5r been  recomputed  uain* 

(1)  A  rectangular  wing  in  steady  flow  at  a  constant 
angle  of  attack  cC  . 
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(2)  A  triangular  wing,  with  leading  edge  sweep 

A.^cat~'a.bV,  in  steady  flow  at  a  constant  angle  of 
attack  ot 

a.  The  Rectangular  Wing 

The  pressure  distribution  on  the  mixed  region  of  this  wing 
is  presented  in  Fig.  VI. 1.  A  comparison  between  these  values 
and  those  from  Table  iV.2  reveals  that  the  inverse-coefficient 
approach  yields  better  accuracy  for  the  pressures  at  the  box 
centers  than  the  direct  method  which  omits  the  effect  of  the 
side-edge  singularity.  However,  the  accuracy  is  poorer  than 
the  direct  method  including  the  refinement  for  the  edge.  The 
discrepancies  between  the  results  of  the  new  method  and  their 
exact  counterparts  can  be  solely  attributed  to  the  approxima¬ 
tion  that  the  pressure  over  each  box  is  constant  and  equal  to 
its  value  at  the  center;  this  is  most  critical  for  a  box  ad¬ 
jacent  to  the  side  edge,  where  the  slope  of  the  pressure  is 
singular .  One  might  be  able  to  get  improved  precision  by  ac¬ 
counting  in  an  approximate  fashion  for  the  known  variation  of 
pressure  over  side-edge  boxes.  Such  a  step  would  require  (as 
in  the  case  of  the  refinement  in  the  direct  case)  developing 
and  tabulating  an  extra  set  of  influence  coefficients. 

b.  The  Triangular  Wing 

The  pressure  distribution  on  this  wing  is  presented  in 
Fig.  VI. 2.  A  comparison  between  these  values  and  the  corre¬ 
sponding  ones  from  Fig.  IV. 13  reveals  that  the  chordwise 
pressure  fluctuations  are  even  more  severe  than  those  en¬ 
countered  with  the  direct  method.  The  dimensionless  total 
lift  and  total  moment  are 


/,=-/.  006, l, 7  ,  A7  =  -/-  £90,99 


These  quantities  are  closer  to  the  exact  values  of  —0-987,c9 
and  -/.3/bj/l  ,  respectively,  than  the  corresponding 

estimates  by  the  direct  method  using  six  chordwise  boxes.  The 
latter  are  -/.o«3, 96  and  -  /.  jys,  /3  from  Table  IV. 4. 

It  has  been  demonstrated  so  far  that  the  inverse  approach 
offers  a  satisfactory  alternative  way  for  estimating  steady- 
state  airload  distributions.  From  the  examples  shown  above 
it  appears  that  one  may  often  expect  to  obtain  better  accuracy 
than  from  the  direct  pressure  influence  coefficients  without 
refinements.  However,  the  new  method  has  a  serious  dis¬ 
advantage  when  extended  to  wings  in  unsteady  motion.  The 
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VALUES  OF 


7 

+/.00000 

/ 

/ 

0.491,51 

0.608,17 

0445/37 

0.391,83 

/ 

0.6/2,19 

0.704,83 

0593,84 

0564.09 

0.393,64 

0.4/0,57 

0326,22 

0.295,17 

/ 

0.671,96 

0  733,25 
0.669,62 
0640,98 

0.527,02 

0.525,55 

0.475,/7 

0.454,37 

0.293,96 

0352,90 

0.248,31 

0.246,75 

/ 

0.703,07 

0.783,65 

0.7/2,07 

0687,43 

0.602,63 

059Q8/ 

0.552,13 

0535,44 

0.403,05 
0.460,54 
0.391,73 
0.39/, 83 

0.290,59 

0.300,0/ 

0233,31 

0.216,35 

/ 

0.734,32 

0805/72 

0. 741,43 
0.7/3,56 

0.653,37 

0634,73 

a  602^5 

0.587,33 

04677/ 

0525,88 

0.470,63 

0.471,02 

0.400,73 

0.40/99 

0.362,09 

0.349,80 

0229,06 

0274,05 

0/91,69 

0/94,96 

Fig.  VI. 1  Dimensionless  Pressures  on  the  Mixed  Region 
of  a  Rectangular  Wing  in  Steady  Supersonic 
Motion  at  a  Constant  Angie  of  Attack.  Com¬ 
parison  Between  Exact  Results  and  Three 
Forms  of  the  Mach  Box  Method 

[Values  given  are  for  the  centers  of  boxes 
in  the  following  order:  (1)  Inverse  pres¬ 
sure  method,  (2)  Direct  method  without  down- 
wash  singularity,  (3)  Direct  method  with 
downwash  singularity,  (4)  Exact.] 
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Fig.  VI. 2 


Pressure  Distributions  on  a  Triangular 
Wing  in  Steady  Supersonic  Flow  at  a 
Constant  Angle  of  Attack  with  the  In¬ 
verse  Pressure  Coefficient  Method. 

( *  * 6 -  cot  A  =  o.W,  M=Jz) 
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form  of  Che  inverse  kernel  K  is  much  more  compllcaced  then 
that  of  iCs  counterpart  F  ,  rendering  the  task  of  evaluating 
the  inverse  pressure  influence  coefficients  a  difficult  one. 
The  downwash  “Ztx.y)  at  a  point  (*,#)  due  to  a  unit  constant 
pressure  on  an  elementary  area  A  inside  the  forward  Mach 
cone  from  can  be  written  in  dimensionless  notation  fcf.  Eq. 
(16b)  of  Ref.  42]: 


V‘  p  [/(WM)3"  (p-y  foty-lfjhjfy 

Wrl 


Eq.  (6.7) 


Compared  to  the  expression  for  the  direct  influence  coefficient 
(Eq.  (2.4)],  the  evaluation  of  the  inverse  coefficient  ac¬ 
cording  to  Eq.  (6.7)  calls  for  an  additional  integration  over 
the  variable  X  ,  Also,  extreme  caution  must  be  exercised  in 
carrying  out  the  >7  -integration,  since  the  integrand  is  highly 
singular  at  .  Although  Eq.  (6.7)  can  be  Integrated  in 

closed  form  when  *= o  to  yield  the  Influence  coefficients 
Eqs.  (6.6a-c),  it  is  a  difficult  task  to  evaluate  (x,y) 

when  A  to  .  Even  for  low  frequencies,  when  K  can  be  approxi¬ 
mated  by  a  finite  power  series  in  &  (cf.  Ref.  42),  one  still 
is  faced  with  the  double  integrations  with  respect  to  ? 
which  require  lengthy  sets  of  computations.  Further  research  is 
recommended  before  any  systematic  tabulation  of  inverse  coef¬ 
ficients  is  undertaken. 

VI .2  The  Lift  Influence  Coefficient 

The  lift  influence  coefficient  is  defined  as  the  total 
lift  over  an  elementary  area  due  to  a  unit  constant  downwash 
over  another  elementary  area.  Consider  two  such  areas  to  be 
the  full  boxes  A,  B  ,  respectively  in  the  Mach  system  as 
shown  in  Fig.  VI. 3. 
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Fig.  VI. 3  The  Use  of  Mach  System  in  Connection  with  the 
Lift  Influence  Coefficient. 

In  order  to  evaluate  the  total  lift  over  A  due  to  B  , 
one  must  first  find  the  pressure  distribution  over  A  .  For 
the  point  a.,  ,  the  influencing  area  is  the  triangle  cde  , 

and  for  the  point  ,  the  quadrilateral  //<V  >€  .  Since  the 

evaluation  of  the  lift  influence  coefficient  requires  a  double 
integration  over  A  ,  it  is  desirable  to  find  the  pressures  at 
all  points  of  A  in  closed  form.  Although  this  is  possible 
for  the  steady-state  case,  the  same  cannot  be  said  when  A  . 
Furthermore  this  method  was  found  to  have  no  particular  ad- 
v*n5age  over  the  ori8*-nal  approach  with  the  possible  exception 
of  smoothing  out"  the  chordwise  and  spanwise  fluctuations  of 
the  airloads  and  allowing  the  specification  of  zero  lift  over 
each  diaphragm  box  rather  than  zero  pressure  at  the  center  of 
each  such  box.  Since  the  usefulness  of  this  approach  is  limi¬ 
ted  to  the  steady- state  case,  and  even  there  the  additional 
complications  far  outweigh  its  slight  advantages,  this  ap¬ 
proach  was  not  pursued  further. 
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APPENDIX  A 

EVALUATION  OF  THE  FUNCTIONS  I„  AND^ 

It  was  pointed  out  in  Section  IV  that  two  additional  func¬ 
tions  roust  be  tabulated  in  connection  with  the  procedure  for 
including  the  effect  of,  the  downwash  singularity  near  a  side 
edge.  These  are  the  I„  -  and  -functions,  defined  by  Eqs. 
(4.17a-b).  The  integrations  required  must  be  carried  out 
numerically,  and  the  recommended  Integration  technique  is  the 
modified  Gauss  quadrature  (Ref.  43),  which  appears  to  give 
satisfactory  accuracy  with  a  minimum  of  computational  labor. 

A.l  Numerical  Integration  Formula 

Since  the  integrand  in  is  a  function  of  z*  only,  it 
represents  an  integral  of  the  type 

~  ■  Eq.  (A.l) 

Applying  the  general  technique  of  Ref.  43  to  Eq.  (A.l),  it 
can  be  shown  that  the  numerical  integration  formula 


Eq.  (A. 2) 


is  exact,  provided  fly)  is  a  polynomial  of  order  or  less 

in  g.  and  the  stations  fa  are  properly  chosen.  The  Hjs  are 
associated  weighting  factors.  The  valu~T  of  ^  are  the  roots 
of  the  equation 


Eq.  (A. 3) 


For  example,  if  /V*  5  ,  these  roots  turn  out  to  be 


fa  ^  O.ozz, 

fa  =  0J87,  831,57V 

fa  =  O.V(,l,S$7t3W 


fa  =  0.748,  33V,  (>60 
0-97S,  V-94  9/o  . 
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The  weighting  factors  Hj  can  then  be  obtained  from  the  fol 


9,H, 


i-  hh  *-  +  hn  =  / 

*  kH.  '  j/,H,  ..  v,Vv-i 

*  ^  ^  **,% +-■  =  ± 


I 

rt-t . ,  N-*,,  Af-/,  ,  N-i  » 

**  ^  *A  H,-f, 


Eqs.  (A. 4) 


o.  t¥9,1HSfj3£t 
Ms  ~  0. 066,67/,  358 


For  N  -  *>  f  one  has 

A/x  =  O.  395,  SZ<J,  Z/s 

Hz  -  0.  269,266,  739 
H3  =  O.  2/9,086,398 

„  den««  the  Integrand  of  the  -Integra] 

If  *  finite  power  series  of  order  (Zfhl)  In  f  ‘  which  takes 
the  same  values  as  /%(**■)  at  the  points  y  -  **  a 

&r&n%w 

~  e  ft  Eq.  (A. 4) 

where  6,  o~  are  those  defined  following  Eqs .  (4.l7a-b). 

s.:-^Sa^S-£“‘4r 

s*r  $i&~B  & S5- Tdsi 
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critical  case  will  be  when  the  frequency  is  maximum.  Let  the 
maximum  reduced  frequency  at  which  this  method  is  to  be  ap¬ 
plied  be  denoted  by  .  In  the  process  of  calculating  the 

aerodynamic  influence  coefficients  for  a  wing  at  ,  one 

need  consider  only  the  €  -range  for  Tn  [Eq.  ( 4.17a)] 

o  c  e  ^  =  <£ _ „ 

For  different  boxes  in  the  diaphragm  region,  both  £  and  cr 
vary.  However,  it  can  be  easily  inferred  from  the  locations 
of  diaphragm  boxes  in  Fig.  IV. 5  that 

<r  ^  ~~  -/ 

€ 

Therefore  for  a  given  e  in  the  tables,  the  only  entries  re¬ 
quired  are  those  for  which 


If  the  Mach  box  system  is  used,  and  at  most  nine  boxes  are 
taken  along  the  side  edge  (beyond  this  number  of  boxes,  the 
inclusion  of  the  downwash  singularity  is  probably  unnecessary, 
as  stated  in  Section  IV),  this  limitation  yields  cr^7.*  There¬ 
fore  the  tabulations  would  have  to  cover  the  ranges 

o  *  <r  <  ,  o  <  7  . 

With  these  ranges  established,  it  has  been  found  that  a 
five-point  integration  will  suffice  for  all  c's  and  cr's  for 
which  g  ■*->  £. 

When  evaluating  the  function  j?  ,  the  situation  is  some¬ 
what  different  because  the  limits  of  integration  are  sfiF  and 
1  rather  than  O  and  l  .  A  numerical  integration  formula 
similar  to  Eq.  (A. 2)  can  be  devised  when  the  lower  limit  is 
some  constant  other  than  zero.  However,  the  formula  will  be 
dependent  on  this  lower  limit.  Since  one  must  tabulate  the 
coefficients  j p  for  several  values  of  <7  ,  an  equal  number 

of  integration  formulas  would  be  necessary.  To  by-pass  this 
complication,  the  following  alternative  is  suggested.  Let 

£„(Zz)  denote  the  integrand  of  the  -integral  (Eq.  4.17b) 


*It  should  be  remembered  that  0~  can  have  only  discrete  values 

when  a  box  system  is  used;  e.g.,  for  nine  Mach  boxes,  (T- ^3,5,7 

JL  3  ...7  j_  s  ...  .  ' 

2 ' 2.  1  2 ' 3  '  3  * 
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£=  e 


(z‘)  Jz 


=  e 


\fr 

Gnfe  )dz 


Eq.  (A, 5) 

For  <  (r  ,  the  arguments  of  the  Bessel  functions  in  Eq. 
(4.17b)  become  imaginary.  Instead  of  taking  the  true  values 
of  Gn  for  2,2<<r  ,  consider 

£,<*').  )i(  );Z( 

-£*’e'3‘"Y  i}>  ***** 

Eq,  (A. 6) 

The  function  G„fez)  is  continuous  across  ,  and  has  also 

at  least  a  continuous  first  derivative.  Therefore  G„(&z)  is 
sufficiently  smooth  to  be  approximated  by  a  finite  polynomial 
in  .  Equation  (A. 5)  then  reduces  to  (since  Gje2)^  Gntez)f>r 

P,  - 


=  e‘><r  {[’gJP)^  -ip.fi!ne  Ue<r*‘dzj 

Eq.  (A. 6) 


Now  both  Integrals  have  the  limits  0  and  /  ;  therefore,  apply¬ 
ing  Eq.  (A. 2)  to  Eq.  (A. 6)  yields 


£- 


,  Leer 


Eq.  (A. 7) 
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The  number  of  integration  points  H  that  i,  i  >  .  be  taken  tr 
ensure  adequate  numerical  integration  depends  oi  he  values  oi 

£  and  cr  .  By  a  method  similar  to  that  for  JL  .  can 
prove  that  tables  of  the  have  to  cover  the  ran  <5.® 

"  *  -e~  ^  '  o  ~  cr-  jg  • 

Once  more^nine  >e  maximum  number  of  boxes  to  b<  taken  along 

the  side  edge.  ^  for  this  case  is 

-  »  A 

X  ^  2. 

It  has  been  found  that  for  the  necessary  tables,  a  .  ve-pcint 
integration  for  £  £  1.6,  and  a  nine-point  integration  for 
1.0  £  £.  3.2  will  suffice,  for  all  ar's  ,  provide*  3,2. 
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